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LASER VELOCIMETER SYSTEMS ANALYSIS APPLIED TO A 
FLOW SURVEY ABOVE A STALLED WING 

Warren H. Young, Jr., * Janies F. Meyers, 
and Timothy E. Hepner* 

Langley Research Center 

SUMMARY 

A laser velocimeter operating in the backscatter mode was used to survey the 
flow above a stalled wing installed in the Langley high-speed 7- by 10-foot tunnel. 
Polarization was used to separate the two orthogonal velocity components of the fringe- 
type laser velocimeter, and digital counters were used for data processing. The veloc- 
ities of the kerosene seed particles were measured with less than 2 -percent uncertainty. 
The particle velocity measurements were collected into histograms. 

The flow -field survey was carried out above an aspect- ratio- 8 stalled wing with an 
NACA 0012 section. The angle of attack was 19.5°, the Mach number was 0.49, and the 
Reynolds number was 1.4 x 10 . The flow field was characterized by the periodic shed- 
ding of discrete vortices from near the crest of the airfoil. 

INTRODUCTION 

An extensive investigation into the process of the dynamic stall of an airfoil at 
large angles of attack, summarized in reference 1, has delineated the features of the 
stalling process. The physical processes involved in recovery from stall and the effects 
of compressibility are considerably less understood. The present experiment was 
designed to determine the applicability of a laser velocimeter to velocity measurements 
in a separated region and to provide a preliminary definition of the flow field about a 
fully stalled airfoil at a Mach number typical of the retreating blade of a helicopter rotor. 

References 2 and 3 describe the use of both hot wires and pitot probes to measure 
the flow in the separated regions on an airfoil at angles larger than the stall angle. The 
limitations on the use of a probe in a separated flow include directional ambiguity, uncer- 
tainty in probe alinement with flow direction, and probe interference effects. However, 
a laser velocimeter offers no flow interference or flow alinement difficulties and is cap- 
able of being modified to eliminate directional ambiguity. The laser velocimeter can 
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also make virtually instantaneous measurements and features a small common volume. 
However, the assessment of the applicability of the laser velocimeter to the measure- 
ments above the stalled wing requires the consideration of the entire laser velocimeter, 
seed generation, and data-acquisition system. 

A fringe -type laser velocimeter was installed in the Langley high-speed 7- by 
10 -foot tunnel. The laser velocimeter used polarization separation to measure two 
components of velocity. Operation in the backscatter mode allowed the measurement 
location (the common volume) to be readily scanned. The scan and velocity measure- 
ments were carried out in a vertical plane through the tunnel center line near a wing 
with an NACA 0012 cross section that was set to an angle of attack of 19.5°. The Mach 
number was 0.49, and the nominal value of Reynolds number was 1.4 x 10 . Because 
the wing was straight, unyawed, and had an aspect ratio of 8, and because all the measure- 
ments were made in the plane of the center span, the flow was treated as two-dimensional. 
Frequencies and time histories were obtained from absolute pressure transducers on the 
wing surface. In addition, interpretation of the laser velocimenter results was assisted 
by observation of tufts in the separated region to determine the local flow direction. The 
wing was completely stalled, and the flow was separated over the portion of the upper sur- 
face that was aft of the crest. 

The use of kerosene smoke for seeding the flow is described in appendix A by 
William V. Feller. An explanation of the significance and interpretation of the statistical 
nature of histograms is given in appendix B by Janet W. Campbell. 

SYMBOLS 

C chordwise position at which normal meets surface, percent chord 

number of velocity measurements in ith histogram interval as percent of D 

D total number of velocity measurements in one histogram 

D diameter of common volume, m 

cv 

E excess (or kurtosis) (see eq. (B4)) 

f signal frequency received from laser velocimeter, Hz 

H spanwise distance measured from the center span, percent chord 
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number of histogram intervals 
histogram interval number 
length of common volume, m 
fringe spacing, m 
free -stream Mach number 

distance from airfoil measured along normal to surface, percent chord 

probability density function 

number of cycles in signal burst 

airfoil Reynolds number based on chord 

skewness of histogram (see eq. (B3)) 

distance along airfoil surface measured from leading edge, m 

period, sec 

trailing edge 

transducer number 

time, sec 

velocities in free -stream and vertical direction, m/ sec 

velocities in coordinate direction 45° below free-stream direction, m/sec 

velocities U/V^ and v/v^ 

velocities in coordinate direction 45° above free-stream direction, m/sec 
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Vp speed of pressure pulses, m/sec 

Vqo free-stream velocity, m/sec 

X,Y coordinate axes (fig. 16) 

x distance upstream from airfoil trailing edge, m 

y vertical height above airfoil trailing edge, m 

a s angle between N and vertical, deg 

a w wing angle of attack, deg 

AU width of histogram interval, m/ sec 

AZ distance between beam focus and common volume, m 

e error, percent 

8 angle between crossing laser beams, deg 

X laser radiation wavelength, m 

a standard deviation, m/sec (see eq. (B2)) 

£2 pressure pulse passage frequency, Hz 

co frequency, rad/ sec 

Subscripts: 

a average or mean as defined by equation (Bl) 

m median 

u direction 45° below horizontal 
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v direction 45 above horizontal 

1 mode or velocity of highest peak in histogram 

2 velocity of second highest peak in histogram 

APPARATUS AND CALIBRATION 
Laser Velocimeter Optics System 

A fringe-type laser velocimeter optics system operating in the backscatter mode 
was used for the tests discussed in this report. The system is described in detail in 
reference 4. The primary advantage of the backscatter mode of operation is the ease 
of maintaining optical alinement between the primary transmitting and collecting optics 
since common optical components serve both functions. Aside from the practical neces- 
sity of operating in this mode in large facilities, the backscatter mode also permits the 
use of simplified optical and mechanical scanning mechanisms necessary for obtaining 
a three-dimensional spatial survey of the velocity flow field of interest. A two-component 
configuration, using polarization separation, was used to measure the velocity compo- 
nents ±45° about a flow streamline parallel to the test-section longitudinal center line. 

The two velocity measurement directions are called U-^ and V^. The laser velocim- 
eter optics are illustrated in figure 1. 

The laser velocimeter depends on the monochromatic and coherence properties of 
the laser. A laser beam is optically split into two parallel beams which are then focused 
with a lens to the point where the two beams cross. In the common volume, at the cross- 
over point, the light waves from each beam interfere constructively and destructively 
forming a fringe pattern. If the two beams are alined in the horizontal plane, the fringes 
appear as vertical planes extending along the optical axis. The separation between two 
adjacent fringe planes is found from 


L„ = 

fr 


2 sin |- 


where A. is the wavelength of the laser radiation and 6 is the angle between the two 
crossing laser beams. As a lira . -sized particle passes through the fringe pattern, light 
is scattered with an intensity proportional to the incident intensity in the fringe pattern. 
If a portion of this scattered light is collected by a photomultiplier tube, an electrical 
signal is obtained whose oscillation rate is the same as that of the scattered light 
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oscillations generated as the particle passes through the fringe pattern. The frequency 
of this oscillation is related to the particle velocity by 



Since the laser beam is linearly polarized, appropriate polarization techniques were 
used to separate the scattered light from a second set of fringes located at the crossover 
point and orthogonal to the original fringes. This second fringe pattern is obtained by 
focusing three parallel laser beams instead of two beams, with each beam located at an 
apex of an isosceles right triangle. The resulting fringe patterns of the two components 
are shown in figure 2. With the two beams on the hypotenuse of opposite polarization 
and the polarization of the third beam rotated 45°, the scattered light from each set of 
fringes formed in the minor legs of the triangle can be separated by depolarizing optics. 
This technique allows the measurement of two orthogonal velocity components. 

The argon-ion continuous -wave laser and the laser velocimeter optical system were 
located in the tunnel plenum chamber to obtain the maximum optical acceptance solid angle 
of the scattered laser radiation. The laser operated at about 2 W and 514.5 nm. The 
laser velocimeter had a focal length of 1.9 m (sufficient to reach the tunnel center line) 
and a collecting solid angle of 0.0458 sr. The common volume, 0.6 cm long with a diam- 
eter of 0.20 mm, yielded 20 fringes with a fringe spacing of 10.1 pm. The optical system 
was placed on a mechanical two -component traversing mechanism to allow movement of 
the common volume in the x- and y-directions. The x-direction was parallel to the flow, 
and the y-direction was perpendicular to the flow in the vertical direction. The move- 
ment was controlled, and the position was calculated by a programmable calculator. This 
entire system was then placed on a pneumatic isolation platform to protect it against 
vibration. The laser was placed in a sealed case vented to atmospheric pressure to pro- 
tect it from the plenum chamber environment during run conditions. The overall assem- 
bly, including the isolation platform, traversing system, laser, and laser velocimeter 
optical system, is shown in figure 3. 

Laser Velocimeter Electronics System 

The interface between the optics system and electronics system consisted of two 
S-20 response photomultiplier tubes, one for each component, and signal conditioning 
electronics. The signal conditioning electronics were composed of line driving amplifiers 
attached to the photomultipliers and low-pass filter banks to remove high-frequency noise. 

The laser velocimeter electronics system, shown schematically in figure 4, mea- 
sures the frequency of the output signal from the optics system, converts it to velocity, 
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and develops velocity histograms from which the mean velocity and standard deviation of 
velocity fluctuations are determined. The histogram data are then stored on magnetic 
tape for offline data reduction. 

High -speed burst counter. - The high-speed burst counter is a device designed to 
measure the period of a high-frequency burst (1 kHz to 100 MHz) of the type received 
from a laser velocimeter. The instrument outputs data in either period measurements 
(time mode) or inverse period measurements (velocity mode). An idealized burst 
received from the laser velocimeter is illustrated in figure 5(a). The pedestal is removed 
by high-pass filters so that the burst is symmetric at about 0 V (fig. 5(b)). A Schmitt 
trigger is then used to convert the burst into a digital pulse train (fig. 5(c)). The first 
pulse in the train is used to clear the counter circuits; the second is used to arm the 
counters; and the third triggers the counters to begin counting pulses from the 500 -MHz 
reference clock. When the tenth signal pulse occurs, the counters are halted. The 
counter now contains the measurement of the period average of the eight previous signal 
cycles based on a reference clock of 500 MHz yielding a period average measurement 
with a resolution of 2 nsec. However, since the counter has a finite size (9 bits for 
8 cycles), frequencies below 3.9 MHz have periods too long to be counted. When this 
occurs, the counter register is shifted to the right by 1 bit, and this shift effectively 
halves the reference clock frequency. The number of shifts is output as a 4-bit expo- 
nent, n + 6. (The exponent bits are biased by 6, e.g., exponent 7 is 1 shift.) For exam- 
ple, if the number 1011010110 exponent 0101 (726 jq exponent 5) were output, the measured 
period would be 726 x 2 x 10' 9 x 2( 5 " 6 ) or 0.726 jusec for 1 cycle which yields a frequency 
of 1.38 MHz. Thus, the digital word from the high-speed burst counter consists of 10 data 
bits and 4 exponent bits. The conversion equation to signal frequency is 

„ 3.2 X 10 10 

Frequency = -= 

2 X Reading 

Since the high-speed burst counter was placed in the velocity mode for the present test, 
the most significant 8 bits of the time reading are used to address a read-only memory 
which contains the appropriate 8-bit conversion to inverse time. The conversion equa- 
tion to signal frequency in the velocity mode is 

Frequency = 7.5 x 2 n x Reading 

These readout schemes permit resolutions of 0.4 percent of full scale in the time mode 
on all exponent scales and 1.6 percent of full scale in the velocity mode on all exponent 
scales. 
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The digital output section in the high-speed burst counter was designed to interface 
with a minicomputer; thus, the digital word was designed to shift exponent scales when 
appropriate to seek the maximum resolution of the measurement. However, the pulse 
height analyzer was designed to accept only 10 bits; this design required the burst 
counter exponent to be set and held for a given histogram. 

Pulse height analyzer . - A pulse height analyzer (PHA) was an instrument used to 
record the number of occurrences of pulses of given amplitudes. In the PHA used in 
this study, the normal analog front end was replaced with a two-channel digital multi- 
plexer designed to accept two 10 -bit words from the two high-speed burst counters. The 
capacity of the memory of the PHA was 2048 channels, a number sufficient for full digital 
resolution (i.e., a 10-bit word had 2^ or 1024 possible values). (In the velocity mode, 
the counter outputs 8 bits which address only the lowest 256 channels.) When the PHA 
received a digital word from the high-speed burst counter, the word was placed on the 
memory address bus. The addressed memory location was then incremented. After a 
large number of words had been received, a histogram was formed in which each channel 
location contained the number of occurrences in that particular digital address (i.e., out- 
put word from the high-speed burst counter, during the measurement sample). The his- 
togram was related to the velocity probability density when the memory addresses or 
PHA bin numbers were converted to velocity by using 

V = L fr (7.5 x 2 n x Bin number) 

where is the laser velocimeter fringe spacing, and the high-speed burst counter is 

in the velocity mode. The histogram was then stored on digital magnetic tape. 

Seeding 

The particles used to scatter laser radiation were generated by a vaporization- 
condensation smoke generator. Kerosene was used as the smoke material. The 
details of the system installation and operation of the smoke generator are found in 
appendix A. 

Since the laser velocimeter measures particle velocity instead of gas velocity, the 
response of the particle to variations in a flow field must be known. The size distribution 
of the particles used to obtain the velocity measurement must be known to calculate this 
response. The distribution was measured in still air with an electro-optical particle size 
analyzer. The instrument measures the size of particles from 0.5 to 8.0 pm. Its prin- 
ciple of operation is based on the light scattering criteria from small particles developed 
by Mie (ref. 5). As the smoke is pulled through the optical system with a small air pump, 
the particles are illuminated with white light, and the scattered radiation is collected over 


8 


a large collecting solid angle at 90° with respect to the illuminating radiation. The 
combination of white light and the large collecting solid angle smooths the Mie intensity 
function (ref. 5) to a straight line on a log-log scale. The collected scattered radiation 
impinges on a photomultiplier tube which yields an electrical pulse of a height propor- 
tional to the particle size. The pulse is then shaped into a rectangular pulse, whose 
height is equal to the peak of the photomultiplier output pulse, and is transmitted to a 
PHA in which the particle size distribution is obtained. A diagram of the system is shown 
in figure 6. The particle size distribution for the generator used in the present test is 
discussed in appendix A. 

Sufficient smoke of a size large enough to be detected by the laser velocimeter 
(above about 3 pm) was generated to yield from 2 to 50 laser velocimeter measure- 
ments per second. The data -acquisition rate depended not only on the smoke -generator 
operation, but also on the part of the flow field under investigation. The penetration of 
the smoke into the separated region of the flow field was noticeably poor. 

Model and Pressure Transducers 

The model wing was installed horizontally near the center line of the Langley high- 
speed 7- by 10 -foot tunnel. The test section was slotted and was about 2 m high by 3 m 
wide. The stagnation pressure was atmospheric, but the stagnation temperature rose 
as much as 30° C above atmospheric because of the fan power dissipation in the closed 
circuit tunnel. 

The model wing is shown in figure 7. The model was an untwisted, unswept, 
untapered wing with an aspect ratio of 8 and 0.1524-m chord. The wingtips were bowed 
down a few mm so that the top center span surface of the wing was accessible to laser 
beams from the windows in the side of the tunnel. The wing section was an NACA 0012 
airfoil, and the surface was smooth aluminum. The wing was mounted so that aerody- 
namic response or buffet was negligible. The wing mount had a natural frequency in 
excess of 40 Hz, and the sting was reinforced. Both the wing and sting were prestressed 
so that no vibration was detectable by means of closeup examination by a motion-picture 
camera. Consequently, it was assumed that all the unsteadiness in the flow was due to 
aerodynamic sources. 

The wing position was slightly upstream of the window which was about 18 cm 
square. The laser beams were therefore angled 15° upstream to reach the wing. The 
field of view of the laser beams was thus cut off from the airfoil surface forward of the 
airfoil crest. Also, when the common volume was behind but near the trailing edge, one 
of the beams intersected the trailing edge after passing through the common volume. 

This condition allowed the location of the common volume to be precisely determined 
when it intersected the trailing edge. 
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The center insert of the wing is shown in figure 8. The upper surface was tapped 
for pressure in nine locations with a 0.7-mm drill. Absolute pressure transducers with 
a silicon diaphragm 3.2 mm in diameter were mounted inside the wing. The transducers 
nearest the trailing edge were mounted with the diaphragm flush with the surface. The 
nominal natural frequency of the transducers was 125 kHz. The drilled holes that trans- 
mitted pressure to the transducers had resonant frequencies above 14 kHz. Within the 
frequency range of this test, the response of the transducers mounted inside the wing 
was found to be as large in magnitude as the flush -mounted transducers. The trans- 
ducers at the 0.978-, 3.955-, 8.927-, 14.925-, 23.952-, 34.988-, 59.985-, 74.905-, and 
95.8-percent chord locations were operational throughout the test. The locations of 
all the pressure transducers are shown in figure 9. 

The pressure transducers were calibrated by installing the entire wing in a vacuum 
chamber. Recordings of transducer output were made for five pressures from 1 to about 
1/2 and for temperatures of about 6°, 17°, 26°, 38°, and 43° C. The transducer outputs 
and a reference voltage were filtered and amplified by dc amplifiers and were recorded 
on FM analog tape. 

The letters in figure 9 denote the location of seven tufts that were mounted on four 
small rods protruding from the airfoil surface. The tufts were mounted at staggered 
spanwise locations from 8-percent to 18-percent semispan from the center line. The 
tufts were not run simultaneously with the laser velocimeter. 

DATA ACQUISITION AND REDUCTION 
Pressure Measurement 

The outputs from the pressure transducers in the model were recorded on 10 chan- 
nels of an FM tape recorder. The signals were amplified and low -pass filtered with 
1-kHz filters to eliminate higher frequencies. Recordings were made at 19 cm/sec 
(7I in/sec) on an intermediate band carrier. The tapes were again low-pass filtered 
at 2.5 kHz and digitized at 5000 samples per second. Because of drift in the zero point 
in the recording system, no valid measurement of the average pressure was possible. 

The digitized records were used to find the standard deviation and the amplitude spectra. 

In addition, the pressure signals were analyzed online by a 100 -bin crosscorrelator. 
These data were taken as photographs of an oscilloscope screen. The delay time between 
bins was 20 psec, except for the autocorrelation of transducer 20 which was taken with 
a 5-psec delay time. The horizontal scope scale was 0.2-psec delay time per division. 
The vertical scale was linear. 
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Laser Velocimeter Data Processing 

The number of reference clock pulses for 8 cycles of the laser velocimeter signal 
from two (Ul and Vl) high-speed burst counters was input to the digital multiplexer 
in the PHA. This allowed two histograms, one each for U L and V L , to be simul- 
taneously determined, placed, and displayed side by side on the 2048-channel (bin) 
a nal yzer — that is, one histogram was contained in bins 1 to 1024 and the second was 
contained in bins 1025 to 2048. Depending on the seeding rate, the time to position the 
laser and to acquire a histogram was about 3 to 10 min. The online determination of 
the velocity interval with the greatest number of measurement occurrences (mode) 
required about another minute. 

The dual histogram data were then stored on a digital magnetic tape. Online, the 
component velocities were estimated by conversion to velocity of the mode. Offline, the 
histogram data were retrieved from the magnetic tape and transferred to the program- 
mable calculator for data reduction and permanent storage on magnetic disk. 

Two pairs of histograms are displayed in figure 10. The horizontal axis is the 
velocity value associated with the velocity intervals. The height of each line is the per- 
cent of measurements that fell into that velocity interval. As explained in appendix B, 
the histogram is an approximation to a probability density function. Appendix B also 
defines the relationship between the quantities determined from the histogram and time- 
averaged values. A synopsis of the physical significance of some characteristic histo- 
gram shapes is also given. 

For each pair of histograms, the average velocity magnitudes, U L>a and V La , 
were calculated for each component. The average or mean velocity in the histogram 
(eq. (Bl)) was physically significant because it was also the time-averaged velocity. 

Three other velocities that were calculated are displayed in figure 10. The median 
velocity, U L m or V L m , has its customary definition; that is, one-half of the particles 
have velocities higher than the median, and one -half have velocities lower than or equal 
to the median. The mode 1 velocity, U L 1 or V L j, is the velocity interval with the 
largest percentage of measurements. The mode 2 velocity, ^ or 2’ was ^eter- 
mined visually and is the speed at the second highest peak in the histogram when that 
peak is higher than its neighbors. For most histograms, there is only one peak. The 
differences between the mean, median, and mode 1 velocities are usually much smaller 
than the illustrative histograms in figure 10. Figures 11, 12, and 13 were constructed 
by assuming that both U L and V L were in the positive direction. This assumption 
defines a velocity vector in the prime direction. Careful examination of the three veloc- 
ity maps reveals small differences near the airfoil surface or in regions of highly 
skewed histograms, but for purposes of examining the overall steady flow pattern, the 
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three are essentially the same. Thus, the overall flow picture is not distorted either by 
the fine details of the histogram shape or by variations in the data-reduction process. 

The higher moments of the histogram were also taken. The standard deviation 
(eq. (B2)) is the root mean square of the variation from the mean or average and is a 
measure of the unsteadiness of the flow. The skew (eq. (B3)) was calculated to indicate 
the degree of symmetry about the average. The excess (or kurtosis) was calculated 
(eq. (B4)) to complete the moments. Excess is a measure of the degree of peakedrress 
of the histogram. A few measurements of skew and excess are available for boundary 
layers in references 6, 7, and 8. 

The fidelity with which the shape of the true probability density function is repre- 
sented by a histogram depends on the total number of measurements D in the histogram. 
Figure 14(a) shows a histogram with fewer than 2 points per bin. Even histograms with 
a sufficient number of points contain by chance empty bins and have a ragged profile 
(fig. 14(b)). 


System Measurement Precision 

The overall measurement precision is obtained by determining the accuracies of 
all variables in the system which affect the accuracy of each velocity measurement. The 
differentiation of f = v/Lf r indicates which variables cause inaccuracies in the velocity 
measurement 

y = M 

2sin| 

d y = A, fdfi \ 

2.0 - . 0 . 0 
Vs in g 2 sin ^ tan j J 

Because of the stability of the laser wavelength, the wavelength is not considered 
a source of error 


dV 

V 


df d0 

f " „ . 0 

2 tan 2 


For the angle 0 used in the present test 
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dV~ df d0 
V ~ f ” 9 

The frequency error df/f is affected by errors in the optics system: parallel- 

ism of the fringe pattern, measurement of the crossbeam angle, and position of the 
common volume. The frequency error is also affected by errors in the electronics 
systems: limiter threshold error, reference clock synchronization error, and quantiz- 

ing error. The crossbeam angle error d e/6 is affected by the accuracy of the physical 
measurement of 9 . 

Optics system . - The three sources of inaccuracies due to the optics system are 
measurement of the crossbeam angle, parallelism of the fringe pattern, and the measure- 
ment dead zone. The inaccuracy in the measurement of the crossbeam angles yields 
a fixed bias error in the velocity measurement. These angles were measured by obtain- 
ing a displacement pattern of the three laser beams approximately 6 m from the common 
volume, locating the center of each beam, and calculating the angle based on the beam 
separation and distance from the common volume. The center of each beam was located 
within an estimated uncertainty of ±1.5 mm, which yielded an uncertainty in angle mea- 
surement of ±0.6 percent. 

If the fringe pattern is not parallel, the fringe spacing is not constant over the 
length of the common volume, and a measurement inaccuracy occurs. (See ref. 9.) 

This problem occurs when the focus of each beam does not coincide with the beam cross- 
-over point. The lack of coincidence is caused by the small divergence in the laser beam 
due to the concave reflectors inside the laser. The problem also adds a bias error since 
the fringe spacing, based on the measured crossbeam angle, applies only when the beam 
focal point coincides with the crossover point. 

The measured beam focal point was located 0.055 m from the crossover point at 
a focal length of 1.9 m. The common-volume diameter 0.20 mm was found by magni- 
fying the common volume and determining the number of fringes within the l/e^ inten- 
sity points and multiplying by the fringe spacing 10.1 pm. The broadening caused by 
nonparallel fringes is determined to be 0.46 percent by the equation 


^cv 

AZ 

LmI 


\ XAZ / 


where L cy is the length of the laser velocimeter common volume, AZ is the distance 
between the beam focus and the common volume, and D cy is the diameter of the 
common volume. 
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The bias error caused by the change in fringe spacing is found by a ratio of the 
diameters since the number of fringes within the common volume is a constant 

e bias = cv — = -2.1 percent 
D cv 

where the beam waist diameter * s §i yen by 



Since the common-volume diameter is always larger than the beam waist diameter, the 
bias error always results in a bias toward the lower velocities. 

The measurement dead zone is caused by trajectories through the fringe pattern 
at large angles from the normal. When a particle does not cross the required number of 
fringes (10) for the data processing instrumentation to determine the signal frequency, 
that data point is lost. This loss occurs if the particle passes too far above or below 
the center of the common volume, or when the angle between the particle trajectory and 
the fringe plane is too small, or when a combination of both conditions occurs. Since 
20 fringes are contained in the 0.20 -mm -diameter common volume of the present system, 
the particle must pass within ±0.086 mm of the center of the common volume with an :!R S 
angle no less than 30° to the plane of the fringes. With the fringes oriented ±45° to the 
tunnel center line, flow angles of less than ±15° are required to obtain simultaneous 
measurements in both components. 

Electronic instrumentation . - The sources of measurement error caused by the 
electronics are: limiter threshold, reference clock synchronization, and quantizing. 

The limiter threshold error (ref. 10) is caused by the time jitter in measuring the 
threshold crossing points due to the Gaussian envelope of the signal burst. This error 
is given by 


e 


1 

12 Q 


96 Q 2 


= 1.06 percent 


where Q is the number of cycles (8) in the signal burst used for processing. This error 
is a positive error which biases the measurement toward the higher frequencies. 

The clock -synchronization error is caused by the integer operation of the reference 
clock counting circuits. At the maximum signal frequency of 20 MHz, the error is less 
than 0.5 percent for a reference clock frequency of 500 MHz. This error is a bias error 


14 


toward the higher frequencies. The quantizing error is similar. This error is caused 
by the 8 -bit accuracy of the high-speed burst counter which yields a bias error of 
0.78 ^ e = 1.55 percent toward the higher frequencies. The 8-bit output also sets 
the resolution of the instrument: 


Resolution, 

m/sec 

Range, 

m/sec 

2.56 

160 to 200 

1.28 

80 to 160 

.64 

40 to 80 

.32 

20 to 40 

.16 

10 to 20 


However, for the purpose of plotting the histograms, the ranges were combined so that 
the histogram resolutions, or velocity interval widths, were 


Resolution, 

m/sec 

Range, 

m/sec 

2.56 

20 to 200 

1.28 

10 to 100 

.64 

5 to 50 


Velocity er ror summary . - The summary of system errors is as follows: 


Velocity error source 

Bias error, 
percent 

Random error, 
percent 

Crossbeam angle measurement 

±0.6 

— 

Diverging fringes 

-2.1 

±0.46 

Time jitter 

+1.06 

— 

Clock synchronization 

+0.25 

±0.25 

Quantizing . 

+1.17 

±0.38 

jTotal error 

-H3.98, - 0.22 

±0.65 


The total effective error is obtained by taking the square root of the sum of the 
squares of the total error. Thus, the total effective error is between 0.69 and 
1. 18 percent. 
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Velocity Measurement Accuracy 

In addition to these system errors, errors caused by common-volume position error 
and particle lag error must be considered for overall measurement accuracy. 

Particle lag . - Since the laser velocimeter measures particle velocities and not the 
gas velocity, the final measurement accuracy is dependent on the ability of the particle to 
follow the flow faithfully. The particle dynamic work performed by several researchers, 
notably Walsh (ref. 11), indicates that the smaller the particle, down to the size of the gas 
molecule, the more faithfully it follows a given flow. However, the particle light scatter- 
ing work performed by Mie (ref. 5) indicates that the smaller the particle, the less light 
it scatters. A compromise was reached by using the laser velocimeter simulation pro- 
gram developed by Meyers (ref. 12). This program showed that the optimum particle 
size for the research conducted during this study was between 1pm and 4 pm. Kerosene 
particles were the seeding material for this study. (See appendix A.) Fewer than 2 per- 
cent of the particles were larger than 3 pm in diameter. 

A test was performed to determine experimentally the response of the particles to 
a flow gradient comparable with the expected gradient in the wing study. The test mea- 
sured the velocity along the stagnating streamline of a hemispherical model at Mach 0.55. 
The experimental results and theoretical prediction are compared with the theoretical gas 
velocity (ref. 13) in figure 15. Based on these results, the small uncertainty due to par- 
ticle lag can be neglected in the flow regions accessible to this investigation. However, 
in the unsteady flow field described in the results section, particle lag errors as large 
at 5 percent for 3-pm particles and 10 percent for 5-pm particles were found by applying 
the methods of reference 10. 

Common-volume location . - The common -volume position was determined using two 
rack-and-pinion driven linear potentiometers for the two mechanical scan components. 

The position accuracy of ±1 mm in the x-component (along the tunnel center line) was 
caused by gear backlash in the rack and pinion. The vertical component also suffered 
from gear backlash but had an additional problem of vertical uncertainty which was caused 
by the action of the pneumatic vibration isolation system. The system's height was depen- 
dent upon tunnel pressure and weight shift, which occurred as the optical system was 
traversed in the x-direction. This problem was removed by touching the model with the 
beam crossover point and obtaining a new zero before making a scan. This resulted in 
inaccuracies of ±1 mm in the vertical component far from the wing. Inaccuracies 
decreased to ±0.3 mm near the model surface. The scans were taken along the line 
normal to the surface of the airfoil as shown by the coordinate system in figure 16. 
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Histogram Interpretation 


Velocity measurements . - The velocity measurements cannot be interpreted without 
consideration of the direction ambiguity. The laser velocimeter measured two compon- 
ents, U L and V L . The positive direction of each of the components is shown in fig- 
ure 16. These two components may be combined to give the velocity. However, since 
only the speed, or velocity magnitude, was measured for each separate component, it 
is not known whether U L is actually in the positive or the negative direction. A similar 
ambiguity is present for V L - When both components are set to positive directions, the 
resultant velocity vector is said to be in the prime direction. The prime direction is 
probably the correct direction for most of the points measured. There are three alter- 
natives to the prime direction: the component alone may be taken in the negative 

direction, the U L component alone may be taken in the negative direction, or both may 
be taken in the negative direction. Thus, the resultant velocity vector has a fourfold 
directional ambiguity. 

Any particle that crossed fewer than 10 of the fringes for a component was rejected 
by the data-acquisition system for that component only. Figure 17 illustrates the fringe 
pattern for the U L component. The situation for the V L component is similar but, 
of course, the fringe pattern is rotated 90°. For the component, a particle trajec- 

tory that made an angle with the fringe plane of less than 30° was not measured. Also, 
trajectories that were too far from the center of the common volume to cross 10 fringes 
were not measured. This accounts for the large disparity in the number of particles 
between the U L and V L component that often occurs. The actual experimental situa- 
tion is complicated by the noncircularity of the common-volume cross section, the uncer- 
tainty in the number of fringes (18 to 20), and the oblate spheroidal shape of the common 
volume. The trajectories are assumed to be straight lines, but since the maximum diam- 
eter of the common volume is 0.2 mm, this is not a serious source of error. This total 
dropout zone of velocities for the U L component is shown in figure 18. The radii repre- 
sent velocity vectors. For the 1.28-m/sec velocity interval resolution, no measurements 
are made in the shaded region. Similar diagrams can be drawn for the V L component 
and for the other velocity interval resolutions. 

Velocity bias. - The histograms can be biased toward the higher velocities by the 
velocity bias effect described in references 14, 15, and 16. This effect is caused by the 
tendency of higher velocity air to pass a greater number of seed particles through the 
common volume than lower velocity air. This effect is based on the assumption that the 
seed particles were homogeneously distributed in the streamline tubes passing through 
the common volume. 

Because the smoke arrived in the test section in streamline tubes, the assumption 
is again made that the seed particles were homogeneously distributed when the smoke 
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entered the fringe pattern. The histograms were corrected for velocity bias by dividing 
the number of occurrences of each velocity interval by the value of the velocity. 

Statistical uncertainties . - Each histogram is a discrete approximation of a proba- 
bility density function. The probability that a value of velocity u lies within an inter- 
val u - to u + was approximated. (See appendix B.) Further, the laser velo- 
cimeter measurement was made over a finite period of time, obtained a finite number 
of discrete velocity measurements, and is thus an approximation of a stationary condi- 
tion. Therefore, there is a degree of uncertainty that the calculated mean velocity and 
standard deviation of velocity are the true stationary mean and standard deviation of the 
flow. The measurement uncertainty in the mean value (from ref. 17) is expressed for a 
95 -percent confidence limit as 

Uncertainty in U„ = ± ~^= 
a 


where a is the standard deviation of the histogram, and D is the number of discrete 
velocity measurements in the histogram. That is, with 95-percent confidence, the true 
stationary mean velocity of the flow lies within a radius of the uncertainty about the mea- 
sured mean calculated from the histogram. Similarly, the measurement uncertainty in 
the standard deviation for a 95-percent confidence limit is 


Uncertainty in 



where E is the excess (or kurtosis) of the histogram. 

In this study, the average uncertainty in all the mean velocity magnitudes measured 
is ±3.0 m/sec or approximately ±3 percent of the local velocity. The average uncertainty 
in the flow angle measurement is ±1.5°. The average uncertainty in the standard devia- 
tion measurement is ±2.1 m/sec or approximately ±6 percent of the local standard 
deviation. 


RESULTS 
Histogram Data 

The histograms for all the measured points are presented in figures 19 to 28. The 
figures are arranged in the order of the chordwise location C at which a line from the 
measurement point, normal to the surface, strikes the surface. The histograms within 
each figure are arranged in order of decreasing distance from the surface N to the 
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measurement point. The airfoil contour and arrows proportional in length to the mean 
velocities for the scan are displayed above each set of histograms. Figure 16 illustrates 
the coordinate systems. Many of the histograms for measurement scans above trans- 
ducer positions 3, 5, and 7 show double peaks. Examination of the histograms for scans 
above transducer position 17 or 19 shows the change of histogram range that was needed 
as the measurement location was moved nearer to the surface. Since the histograms are 
considered raw data, no correction for velocity bias has been applied. 

The aerodynamic conditions corresponding to each measurement point are recorded 
in table I. Although the Mach number was close to 0.49 throughout the test, the free- 
stream velocity and Reynolds number varied from day to day. The nominal value of 
Reynolds number was 1.4 x 10 6 and the nominal value of free -stream velocity was 
170 m/sec. The day to day variation in free -stream velocity was less than 10 m/sec. 
Occasionally, because of the fixed upper limit on the velocity of the laser velocimeter 
measurements, a change in tunnel velocity caused some data to be lost. (See fig. 19, 
points 1 and m.) 

The wing angle of attack changed from 19.44° to 19.50° as the cables from the 
sting to the floor stretched (fig. 7). Table I identifies the location in the x,y and in 
the C, N coordinate systems. In addition, the transducer identification number T r 
that corresponds to the chordwise position C is tabulated. The angle a s between a 
normal to the surface and the vertical is also presented. The histogram data corrected 
for velocity bias include the U L and V L components of average velocity, median 
velocity, mode 1 velocity, standard deviation, skew, excess, and number of particles. 

The histogram characteristics are summarized in figures 29 to 38. The consis- 
tency of the measurements is poor for transducer positions 3 and 5 but improves con- 
siderably for the subsequent transducer positions. The consistency of the measurements 
is affected by changes in the velocity range of the laser velocimeter and by changes in 
the free -stream velocity for cases when there were significant numbers of particles 
near the upper limit of the range. The measurements at low average velocities (below 
about 40 m/sec) are affected by the low velocity dropout zone. The scatter is expected 
to increase for the higher moments of the histograms. However, definite trends can 
often be detected in skew and, occasionally, can be detected in excess. The plots of 
standard deviation indicate that there is considerable unsteadiness in the flow. Values 
of 25 to 50 m/sec of standard deviation are seen in most of the plots. The free-stream 
velocity component in both the and directions is Vrj, cos 45° or about 

120 m/sec. The magnitude of the velocity fluctuations is a large fraction of the free- 
stream velocity. 
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Pressure Correlations 


The outputs of the pressure transducers show evidence of large, nonramdom fluctu- 
ations in the flow. The crosscorrelations and autocorrelations of the pressures are shown 
in figures 39 to 45. The pattern of correlation among the transducers indicates a moving 
pressure pattern. The correlations for transducer 9 in figure 42 afford a good example. 
The picture for 9D1 (9D1 indicates the crosscorrelation between transducers 9 and 1 
with the delay on 1) shows a relatively weak correlation that peaks at bin 36. Bin 36 
corresponds to a delay time At of 720 psec. This correlation means that the pressure 
variations at transducer 1 resemble the pressure variations that occur 720 p sec later 
at transducer 9. If the cause of the pressure variations was moving downstream from 
transducer 1 to transducer 9, the speed V p could be estimated by dividing the distance 
between the two transducers by the delay time. 

The pictures for 9D3, 9D5, and 9D7 exhibit a much closer correlation than 9D1. 

The delay times decrease in a systematic manner from transducers 3 to 5 to 7. The 
delay time between transducers 3 and 5, for example, can be obtained by subtracting the 
delay 9D5 from 9D3. A more reliable value is the crosscorrelation 5D3 shown in fig- 
ure 39 or the one shown in figure 40. The delay times are summarized in figure 46. 

The autocorrelation of transducer 9 (labeled 9D9) is shown in figure 42(e). The 
shape is characteristic of a sinusoidally varying pressure with some variation in period. 
The first peak indicates a delay of 240 psec or a frequency of 2083 Hz. If the frequency 
is associated with the passage of the pressure pulse at speed V p , the spacing of the 
pulses can be determined. The results for all transducers are displayed in figure 46. 

Figure 42 shows that for transducers downstream of transducer 9 there is a decrease 
in the degree of correlation. The correlation 19D9 does not allow the calculation of a 
reliable delay time. 

Figure 46 summarizes the correlations of the pressure transducers. The delay 
times in figure 46(a) are the times required for the pressure pulses to move from the 
upstream transducer to the downstream transducer. Since the distance between the two 
transducers is known, the average speed of the pressure pulse over the distance between 
the transducers can be calculated. These speeds are displayed in figure 46(b) for adja- 
cent pairs of transducers. The points are plotted at a surface distance s halfway 
between the two transducers. 

Additional information is available from the autocorrelations. Figure 46(c) shows 
the frequency of the response of each transducer to the pressure pulse passage. The ratio 
of the speed of the pulse V p to the frequency of the pulse passage £2 gives a spacing 
for the pulses. The spacing, shown in figure 46(d), is calculated from the effect of the 
pressure pulses on the surface. Since the pulse speed is calculated halfway between the 
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transducers and since the frequency of passage is known at each transducer, the frequency 
of passage was interpolated to find the pulse spacing. 

INTERPRETATION OF RESULTS 


Vortex Modeling 


The double peaks in the histograms of transducers 5 and 7 are shown in figures 20 
and 21. The discussion of histogram shapes in appendix B indicates that a periodic vari- 
ation of velocity in time is characterized by a two-peaked histogram. 

The correlations of pressure data taken together with the double peaks seen in the 
histograms suggest that a series of vortices is shed periodically from near the crest of 
the air foil and that these vortices move downstream at an increasing speed. The phenom- 
enon of vortex shedding from bluff bodies has been well documented at lower Reynolds 
numbers (ref. 18) and has been named as a cause of noise generation on airfoils at higher 
Reynolds numbers (ref. 19). The present case was examined for the presence of a region 
through which vortices move downstream. A model of a vortex region was constructed 
and was programmed on the calculator. Simulated histograms were then constructed by 
using a Monte Carlo approach with the calculator by simulating the laser velocimeter 
method of histogram generation as closely as possible. 


The vortex model is illustrated in figure 47. The vortex region extended 10. 5 cm 
upstream and downstream of the point at which the histogram was to be calculated even 
though the velocities were dominated by only one or two vortices near the center. The 
vortices moved parallel to the line chosen as the center line. The distance above or 
below the center line was random with a standard deviation of 1 cm. The vortex spacing 
was set at 3 cm in the direction parallel to the vortex path center line which corresponds 


ated a tangential velocity given by 
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where r is the radial dis- 


chosen by trial and error. Although only a qualitative reproduction of the double peak 
histograms was possible, the value of circulation was probably better defined than the 
other empirical values in the vortex model. When compared with an equivalent airfoil 
lift coefficient (about 0.23), the vortex circulation did not seem excessive. However, the 
circulation was sufficient to reproduce the large velocity difference in the two peaks of 
the histograms above transducer positions 5 and 7. 


In the application of the vortex model to the simulation of the histograms in the scan 
above transducer position 7 (fig. 48), the center line of the vortex path was taken to be 
horizontal and was fixed in height at the location y = 0.0673 m which is the height of the 
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measurement point labeled "f" in figure 48 or run 39 in table I. Above the vortex path 
center line, the simulated histograms showed velocities far above the free -stream value. 
The free-stream velocity was 170 m/sec, so the two components of the free stream in 
the and directions were about 120 m/sec. The measured histograms at 
point "a" also showed the high velocities. At point "a," both the measured and simulated 
histograms showed a relatively small standard deviation. This small deviation indicates 
a flow that has only small fluctuations. At points "b," "c," and "d," the measured histo- 
grams showed a dominance of velocities below the free-stream value. For the simulated 
histograms, the below -free -stream velocities did not appear until point "d" was reached. 
This indicated some defect in the geometry of the vortex model. At points "e" and "f," 
both measured and simulated histograms showed double peaks and an extremely large 
standard deviation. The velocities at which the peaks occur do not match exactly, but 
the trend of increasing dominance of the lower peak is clear in both the simulated and 
measured histograms. Points "g" to "j" showed very good correspondence in shape 
between the measured and simulated histograms. At these points below the vortex path 
center line, the double peaks changed to a single low -speed peak below free-stream 
velocity. The measured histograms nearest the surface were not simulated because the 
model did not include viscous forces. 

A closer correspondence between the measured and simulated histograms could 
probably be achieved by empirical adjustment of the geometry and strength of the vortex 
model. This was not attempted primarily because a quantitative simulation would require 
the inclusion of a prohibitive number of factors. The present model does not account for 
the expansion of the flow about the model, viscous effects, free-stream fluctuations, 
differences in the U L and components, curvature of the vortex path center line, 

the starting point of the vortices, random variations in vortex strength (although random 
variations in location were simulated), nor any variation in core size. The simulated 
vortex spacing of 3 cm along the center line was based on the results of the pressure 
crosscorrelation measurements. Vortex spacing was the only parameter in the vortex 
model that was not set by trial and error comparison of the measured and simulated 
histograms. 

The histograms for the scan above transducer position 5 exhibited double peaks 
near the surface. As shown in figure 30 and figure 49, the points away from the surface 
showed a much higher velocity for the Uj_, component than for the V^ component. 

This characteristic was simulated by using exactly the same vortex model except for a 
15° tilt of the center line to the free-stream direction. The pattern of the measured 
histograms for transducer position 5 was not as clear as for transducer position 7. The 
simulation shows enough correspondence, however, to support the supposition that the 
nonsymmetric appearance of the measured histogram was caused by a tilt of the vortex 
path center line. 
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Vortex Locations 


The model vortex path center line is very close to the surface at transducer 5, but 
it is much higher from the surface at transducer 7. This large vertical movement of the 
vortices is not consistent with the time -averaged velocity flow field. However, it is 
possible that at the time of the vortex passage, the instantaneous flow field could have a 
large vertical component. This hypothesis can be illustrated by the oscillating -airfoil 
dynamic stall work in reference 1. Another illustration of the possibility of large vertical 
velocities is found by combining the mode 1 and mode 2 peaks of this investigation. Since 
the histograms lack time information, the mode 2 peak of one component might occur at 
the same time as either the mode 1 or mode 2 peak of the other component. If both histo- 
grams have two peaks, there are four such possible combinations, as shown in figure 50. 
Thus, the instantaneous vortex convection velocity may not resemble the mean velocity. 

Similar simulations of the histograms downstream of transducer position 7 were 
not attempted. Many of the histograms for transducer 9 (fig. 22) were not sufficiently 
well defined. Also, the vortex model assumptions of inviscid vortices with fixed spacing 
became less valid for the downstream flow. However, double peaks were also observed 
in the trailing-edge region. The points labeled "r," "s," and "t" in figure 28 were for 
flow from the lower or pressure side of the airfoil. This indicates that vortices were 
shed from both the suction and pressure sides of the airfoil. The transducer near the 
pressure side trailing edge, transducer 20, was crosscorrelated to transducer 7 and to 
transducer 3 as shown in figure 45. The delay time for a positive correlation was very 
long (between 130 and 160 p sec). Although the cause or effect cannot be discerned, the 
lower and upper surface unsteady pressures are definitely correlated. 

Frequency Analysis 

The photographs of figures 39 to 47 were digitized, and the waveforms were pro- 
cessed by a conventional computer fast Fourier transform. The real part of the Fourier 
transform is presented in figures 51 to 57. The frequency resolution of the transforms 
is 500 Hz only. The thickness of the spectrum peaks would probably be much less if a 
higher frequency resolution were possible. The correlations show very little high- 
frequency content, and there is generally one dominant frequency in the spectrum. 

The Fourier transform of a pressure autocorrelation is the power spectrum of the 
pressure. The narrowness of the frequency peak for the pressure transducers near the 
airfoil crest indicated that the vortex initiation occurs at very regular time intervals. 

The narrow peaks of the spectra for the pressure transducers near midchord imply that 
the convection speed of the vortices varies little from vortex to vortex. The extreme 
regularity in time of the vortex shedding is reminiscent of the shedding of a Karman vortex 
street at much lower Reynolds numbers on a circular cylinder. 
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The power spectrum of the pressures was also found directly from the analog tape 
recording. The results in figure 58 have been analyzed from 0 to 2500 Hz with a band- 
width of 9.77 Hz. Above 1 kHz, a large roll-off (an effect of the low-pass filter) can be 
seen by comparison with figures 51 to 57. Consequently, only transducers 20 and 2 show 
the peaks found in the unfiltered autocorrelation data since the peak for each transducer 
is about 1160 Hz. The frequency on the lower surface is about one -half of the frequency 
on the upper surface. 

The peaks of the Fourier transforms of the cross correlations could be broadened 
by changes in the timing of the vortex passage just as the Fourier transforms of the auto- 
correlations are broadened. In addition, the peaks could be broadened by changes in the 
magnitude of the pressure pulses. Since the strength of vortices is conserved in poten- 
tial flow, the most likely explanation of the broadening of the peaks of the crosscorrelation 
is a variation in the height of the paths of the vortices. From figures 51 to 56 the path 
irregularity along the length of the chord can be seen. The vertical position of the hori- 
zontal axis of some of the scope pictures could have been affected by the instrument mem- 
ory overflow. In addition, the unusually large dc component in several crosscorrelations 
could be caused by this phenomenon. 


Vortex Paths 

The vortices that pass over the upper surface of the airfoil cause large variations 
from the average velocity. Figures 59 and 60 can be used to compare the average veloc- 
ity vector and the two components of standard deviation. Regions of high velocity and 
large standard deviation can be discerned in the flow field. To aid in the visualization of 
these regions, the resultant average velocity vector was converted to Mach number by 
assuming that the sonic speed was frozen at its free -stream value. Figure 61 shows the 
contours of constant Mach number that result from this representation of average velocity. 
The figure was constructed by means of a spline fit on a 1 -mm grid. A tolerance in Mach 
number of ±0.002 was applied to each contour line. 

Two significant regions of high speed can be seen. The highest Mach numbers 
occur above the midchord of the airfoil. This could occur because of an expansion about 
the separated region combined with the action of the vortex sheet. The high velocity 
region just behind the crest of the airfoil was also seen in the vortex simulation for the 
scan above transducer position 7 (fig. 48). 

A similar procedure was followed for drawing Mach contours of the standard devia- 
tion. Figure 62 shows the regions of unsteadiness in the flow. Three significant regions 
of high standard deviation are apparent. One is near the midchord, and one is above the 
crest. Relatively large values of standard deviations should occur near the center line 
of the paths of the vortices shed from near the crest. 
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The center line of the vortex path can be expected to be characterized by relatively 
high values o£ standard deviation of velocity and zero histogram skew. To track the vor- 
tex paths from the two locations found by vortex modeling for transducer positions 5 
and 7, figure 63 was constructed. In addition to positions of relative maxima of the 
resultant standard deviation, the locations of low values of skew are shown. Also plotted 
are the points where the average velocity component along the chord equals the pressure 
pulse velocity found from the crosscorrelations of pressure. The criteria chosen for the 
path center line did not coincide to define a single track. The center line of the path of 
the vortex system might also influence the average velocity values. The location of 
relative minimums of average velocity with respect to y is also shown in figure 63. 
Except for the simulations at transducer positions 5 and 7, the attempts to fix the center 
line of the vortex paths have not been successful. This failure may be due to the exis- 
tence of a multiplicity of vortex tracks, misapplication of the criteria chosen to define the 
center line, or obfuscation of the data by influences such as viscosity or flow expansion. 

Before proceeding to consideration of viscous -dominated histograms, the inferences 
drawn from the vortex modeling are summarized. The following statements were based 
on the correlations of the pressure measurements, the well-known phenomenon of vortex 
shedding from bluff bodies, and the simulation of the histograms by the vortex model: 

(1) A series of discrete vortices are shed from near the crest of the airfoil at 
regular time intervals. 

(2) The vortices initially move upward, but the movement is nearly in the free- 
stream direction by the chordwise position of 15-percent chord. 

(3) The vortex speed accelerates from about one -half of the free -stream velocity 
near the crest to nearly free -stream velocity near the trailing edge. 

(4) The vortex paths become somewhat less regular toward the trailing edge, but 
the repeatability of the time of the vortex passage remains high from leading edge to 
trailing edge. 

(5) A second region of vortex shedding originates on the lower surface of the airfoil. 

Histograms for Small Average Velocities 

A consideration of the overall flow picture as shown by figures 59 and 60 suggested 
that the flow near the surface behind the midchord of the airfoil should be reversed. 

This reversal was confirmed by the behavior of the tufts. The location of the tufts is 
shown in figure 9. The tufts were observed by a television camera closeup and were 
recorded on video tape. For several tufts, a predominant direction could be discerned. 
The flow at locations A, B, and F of figure 9 was definitely in the downstream direction. 
Tufts C, D, E, and G oscillated between the upstream and downstream directions, and 
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tufts C and D were particularly unsteady. The dominant direction of tufts D, E, and G 
was in the reversed flow or upstream direction. Thus, the average velocity, as plotted in 
figure 36, for example, should have been negative in the region from y * 0.03 to 0.05 m. 

The laser velocimeter measurement of absolute velocity components, combined 
with a range that did not extend to zero, precluded the direct detection of negative veloci- 
ties. However, an attempt was made to infer the flow-reversal region and to assess 
the effects of the measurement dropout zone near zero velocity. 

An interpretation of the minimum velocity calculated from the measured histograms 
can be given from figure 64 which illustrates a Gaussian velocity probability density 
function with an average velocity of 30 m/sec. The full probability density function is 
drawn in figure 64(a). The present laser velocimeter measured only the magnitudes of 
velocity so that the part of the curve to the left of the zero velocity is folded over to the 
right and added to the positive velocity probabilities. Figure 64(b) illustrates this effect 
combined with a dropout zone of ±5 m/sec. The 5 m/sec dropout zone in the laser velo- 
cimeter is caused by a high-pass filter used to remove the laser velocimeter signal 
pedestal so that figure 64(c) is the probability density function that the laser velocimeter 
would have constructed. The theoretical probability density function had a mean value 
of 30 m/sec and a standard deviation of 40 m/sec. The mean in figure 64(c) is 44 m/sec, 
and the standard deviation is 24 m/sec. This increase in average velocity and decrease 
in standard deviation were seen in figure 36. 

Figure 65 reproduces the average velocity and standard deviation for transducer 
position 17 histograms. The symbols in figure 65(a) are measured values. A hypothet- 
ical velocity profile, labeled profile 1, has been superimposed on the measured velocities. 
The average measured velocities above 40 m/sec have been assumed to be unaffected by 
the histogram dropout zone, and the hypothetical profile closely coincides with the mea- 
sured values. The corresponding standard deviation profile is also labeled profile 1. 
Ass um ing the velocities followed a Gaussian probability density function permitted use of 
the computer to generate simulated histograms. The values of average velocity and 
standard deviation (the hypothetical profile 1 values) were input to the computer simula- 
tion. The computer was programmed for the dropout zone and velocity interval resolution 
of the measured histogram. Thus, the computer simulated the distortion that the laser 
velocimeter dropout zone and resolution introduced into the measured histograms. The 
average velocity was assumed to be in the free -stream direction. Once the simulated 
histograms were generated, they were used to calculate a simulated average velocity and 
standard deviation. The computer results are shown as symbols in figure 65(b). The 
average velocity reaches a minimum value of about 20 m/sec and levels off. The three 
values of U L at y = 0.0298, 0.0352, and 0.0396 m are smaller than the Vl com- 
ponent because they were taken with a smaller velocity range. The standard deviation 
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exhibits a similar behavior. This similarity shows that the simulated average velocity 
and standard deviation (the symbols in fig. 65(b)) are a reasonable match for the measured 
values (the symbols in fig. 65(a)). Thus, the dropout zone affects the hypothetical veloci- 
ties given by profile 1 by causing the simulation to give the measured results. The 
simulated histograms could be considered a satisfactory approximation (fig. 66) to the 
measured histograms when the uncertainty in the real probability density function, the 
approximation to the exact size and shape of the dropout zone, and the possible bias 
toward higher velocities (ref. 12) have been taken into account. 

The low velocity dropout zone and the measurement of velocity magnitude apparently 
destroy the possibility of deducing average velocities below 40 m/sec. To confirm this, 
a second profile for transducer position 17 data was drawn. Figure 67(a) illustrates the 
assumed profile, and figure 67(b) shows the resulting average velocities and standard 
deviations. The changes from profile 1 to profile 2 are nearly insignificant. Figure 66 
illustrates the measured, profile 1, and profile 2 histograms. No satisfactory basis 
could be found for preferring one set of simulated histograms over the other. Therefore, 
it was concluded that any further analysis of the histograms with an average velocity 
below 40 m /sec was futile, and that histograms with larger average velocities were little 
affected by the low velocity dropout zone. 

EVALUATION OF SYSTEM PERFORMANCE 

The experimental system comprised the laser velocimeter and data-acquisition 
system, the model and model instrumentation, the wind tunnel, and the seeding. The 
advantages of the system are readily apparent and need little analysis. The major defects 
are peculiar to the present investigation and are correctable. 

The laser velocimeter optics system yielded velocity measurements over a small 
volume with no flow interference. These features were vital for the successful measure- 
ment of the separated flow above the model. The principal disadvantages of the laser 
velocimeter were the dropout zone that affected the mean velocity below 40 m/sec and 
the minimum seed size (about 3 pm) visible to the system. A Bragg_ cell and increased 
optical efficiency will correct both defects. The scan system was satisfactory but was 
limited by the size of the window and the fixed viewing angle of the model. 

The data-acquisition process was never limited by the accuracy of the measurement 
of the seed particle velocities. Even the principal inaccuracies, which were introduced 
by interpretation of the individual velocity measurements, had an inconsequential influence 
on the conclusions of this investigation. The histogram velocity bias, the lag between gas 
and particle velocity, and the statistical uncertainty in the histogram moments are the 
areas presenting the greatest opportunity for improvements in accuracy. The statistical 
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uncertainty can be better controlled, and the velocity bias can be eliminated by an auto- 
mated data-acquisition system that records the value and time of each measurement. 

The time information will also allow frequency spectrum analysis of the velocities. The 
particle lag is calculable for flows for which the time history is known by the methods 
in reference 11. 

The data-acquisition rate was controlled by the seeding of the flow. Improvements 
in both data-acquisition rate and particle lag minimization could be achieved by control 
of the smoke injection position and smoke particle size. Injection of the smoke directly 
into the separated region would be desirable. 

The instrumented model had two major deficiencies. Definition of the magnitude 
of the pressure measurements would be valuable. A larger model chord would have the 
advantages of more realistic Reynolds number simulation and reduced error (relative to 
model size) because of the control volume position uncertainty and control volume size. 
The point of closest approach to model was limited by the diffuse reflection from the 
wing surface. The effects of diffuse reflection can be minimized by smooth model sur- 
faces and improved error rejection in the signal processing electronics. 

CONCLUSIONS 

A survey of the flow field above a stalled wing was conducted using a laser velocim- 
eter that used an argon-ion laser, a backscatter mode of operation, and polarization sepa- 
ration of the two components of velocity. The data acquisition, including the interpreta- 
tion of the histograms of velocity, was analyzed for accuracy. The model was an 
aspect -ratio -8 wing installed in the Langley high-speed 7- by 10 -foot tunnel at an angle 
of attack of 19.5°, a nominal Reynolds number of 1.4 x 10^, and a Mach number of 0.49. 
Unsteady pressure measurements were made on the wing surface. An analysis of the 
results has led to the following conclusions: 

1. The laser velocimeter has measured the components of high velocity smoke 
particles with less than 2 -percent uncertainty. The accuracy of the mean velocity data 
was limited by the data-acquisition and data -reduction process rather than by the mea- 
surement process. 

2. Velocity measurements below 40 m/sec cannot be made by the present technique. 
The removal of the dropout zone, and incidentally the fourfold ambiguity in direction, by 

a Bragg cell is the preferred solution to this problem. 

3. A flow pattern has been deduced which is consistent with the pressure measure- 
ments and the histograms. It does not explain all the features of the data and must be 
regarded as tentative. The primary feature of this flow pattern is a series of discrete 
vortices that are shed from near the crest of the airfoil at regular time intervals. The 


28 



vortices initially move upward away from the airfoil surface. The movement is very 
nearly in the free -stream direction by the chordwise position of 15-percent chord. The 
vortex motion accelerates from about one-half the free-stream speed near the crest to 
near free-stream speed near the trailing edge. The vortex passage time occurs at very 
regular intervals although the trajectory heights are somewhat irregular from vortex to 
vortex. A second set of shed vortices originates from the airfoil lower surface. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 10, 1977 
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SEEDING 

William V. Feller 
Langley Research Center 

Kerosene was chosen as the flow seeding material as a compromise between the 
requirements of low vapor pressure at tunnel conditions to reduce particle size loss by 
evaporation, and high enough vapor pressure to avoid accumulation of liquid on the tunnel 
walls. Because kerosene is flammable, the potential fire hazard had to be evaluated. 

The concentration of kerosene vapor in air at its lower explosive limit is 0.7 per- 
cent by weight (ref. 20). One-fifth of this value was chosen as the maximum concentra- 
tion allowed to accumulate in the tunnel. When the air exchange during the tunnel cycle 
is neglected (a conservative approach) and when evaporation of all kerosene added for 
seeding is assumed, calculation from the tunnel volume (8000 m3) and mass of air 
(9750 kg) indicates that the maximum total amount of kerosene that can be introduced 
is 13.6 kg. At an estimated input rate of 3.78 liters (3.17 kg) per hour, seeding for a 
period of 4.3 hours would be permitted before reaching the limiting concentration 
selected. After this time, the kerosene would be shut off, and the tunnel would be run for 
5 minutes with maximum air exchange to reduce the concentration and to evaporate any 
local accumulation of liquid kerosene. The calculation outlined is conservative and 
neglects the continuous air exchange during flow. Consequently, the fire hazard is 
believed to be acceptably low. 

The particle production and distribution system was mounted on the tunnel turning 
vanes 25 m upstream of the test section. (See fig. Al.) Earlier smoke visualization tests 
had demonstrated the need for distributing the smoke over an area because the streamline 
tubes whip around in the test section with an amplitude greater than the width of the 
smoke plume from a single ("point") source, giving an intermittent laser velocimeter 
signal. 

The smoke distribution apparatus shown in figure Al consists simply of a half- 
cylinder 3.8 cm in diameter upstream of each outlet tube. The half-cylinder extends 
laterally about 1 m. The cylinders are fastened to the turning vanes and are, therefore, 
swept 45° to the flow direction. The smoke produced at the vapor outlet tube is trapped 
in the wake flow of the swept cylinder and is spread laterally as the smoke escapes 
through the wake-stream mixing region. A flow rate of about 0.3 kg per hour was found 
to be sufficient to generate 2 to 50 laser velocimeter measurements per second. 

The smoke (particle) generator (fig. A2) was a Howerton evaporating-condensing 
system that was mounted on the tunnel turning vanes. The kerosene is driven from a 
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pressurized reservoir through three parallel tubes 3 m long (1.8-mm o.d., 1.3-mm i.d.) 
wrapped helically around a 600 -W electrical cartridge heater. The kerosene is vaporized 
in the tubes. On mixing with the tunnel airflow, the kerosene vapor condenses to droplets. 
Size distributions taken before tunnel installation and measured by an electro-optical par- 
ticle size analyzer are shown in figure A3. (The distribution is normalized to the peak 
value and is cut off below 0.5 pm because of limitations of the particle size analyzer.) 

The median particle size is below 1 pm and fewer than 2 percent are larger than 3 pm. 
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Figure Al.- Installation of smoke generator in Langley high-speed 7- by 10-foot tunnel. 
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Figure A2.~ Disassembled smoke generator. 
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APPENDIX B 


INTERPRETATION OF HISTOGRAMS 

Janet W. Campbell 
Langley Research Center 


At any fixed point in the airflow, an instantaneous velocity measurement at time t 
is considered a sample of the time-varying velocity U(t). However, neither the exact 
measured value of U(t) nor the value of time at which the measurement was made is 
recorded. Instead, the fact that the measured velocity lies in the velocity interval 
Uj ± is recorded by incrementing the count in a bin labeled "i." No record is kept 
of the exact value U(t), and the time t is thereafter unknown. The number of measure- 
ments D is large so that the length of time required to make these measurements is 
much longer than any of the periodic variations of the velocity. An assumption is made 
that the velocity U(t) is repetitive in time in the sense that the percentage of measure- 

W 

ments in each bin is stationary in time. The number of counts in bin i as a percentage 
of total number of measurements is denoted as Cp The distribution of Ci is thus 
assumed to be stationary or constant in time. This assumption allows the measurements 
to be analyzed by the statistics of a random variable and allows the values of Ci to be 
assembled into an approximation of a probability density function that is called a 
histogram. 

A probability density function P(U) is the mathematical tool used to define the 
probability that a value of U lies in an interval U - ^ to U + It is by definition 
a function that satisfies the following criteria: 

P(U) ^ 0 



1 


and 



dU = Probability that 


U is between 


and 


U 


k 


The histograms constructed from the experimental measurements are approximations 
to a probability density function. To construct each histogram, the range of observed 
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velocities was divided into I intervals (or bins) of equal width AU with midpoints 
Uj = (i - AU. The fraction of the velocities falling into the ith interval, to 

was counted and denoted Cp Then the probability of the velocity lying between 

Ci 

Uj - and U| + is estimated by Cp where 

Ci Ct 

l c i - 1 

i=l 

The value of Cj is an approximation of P(Ui)AU where P(Ui) is the probability 
density evaluated at the center of the ith velocity interval, and AU is the width of the 
velocity intervals used to construct the histogram. 

In addition to the histogram, other quantities used to describe the probability den- 
sity of the velocity are its mean (or average) standard deviation a, variance o^, 
skewness Sr, and excess (or kurtosis) E. The mean of a velocity probability density 
is simply the average velocity defined by 

U a = T UP(U) dU 


The mean is estimated from a histogram by 

I 

U a - ^ CiUi 


i=l 

The mean of the histogram can be shown to be equal to the time averaged velocity 
which is given by 

T 

U = lim ^ f U(t) dt 

X— oo 1 ^0 


(Bl) 


U 


The equivalence of the histogram mean and the time average can be visualized by 
considering the construction of a histogram from a time -varying signal. In sketch (a), 
the spaces between horizontal lines represent the velocity intervals or bins. The dis- 
tance between the horizontal lines is the velocity interval or resolution AU. The line 
representing U(t) has been dotted to show points at which measurements were made. 
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Velocity 



Sketch (a) 


In this example, a total of 47 measurements was made. The bin at the bottom, which 
corresponds to the smallest velocity, contains a count of four measurements. The Cq 
for that velocity interval is therefore 0.085 or 8.5 percent. The central or third bin has 
the largest count. The velocity associated with this bin (2.5 AU) is called the mode 
velocity. To find the mean velocity, the velocities Up weighted by the values of Cp 
are averaged 


U a = 0.085U! + 0.170U 2 + 0.340U g + 0.255U 4 + 0.148U 5 = 2.706AU 

The mean velocity of the histogram does not necessarily occur at the highest peak or 
mode. The mean reflects the amount of time spent by U(t) in each of the bins or 
velocity intervals. This is what makes the histogram mean or average equal to the 
velocity average in the time domain. 

Further information on the time dependence of U(t) can be extracted from the 
histogram. The variance of a probability density function is the average squared distance 
from the mean. It is defined by 

a 2 = f (U - U a ) 2 P(U) dU 

The variance is approximated from a histogram by 

I 

<r 2 = J Ci(Ui - U a ) 2 (B2) 

i=l 
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Since the variance has units of velocity squared, a more physically meaningful measure 
of the dispersion of the velocity distribution is a, the square root of the variance a 2 . 
The quantity a is called the standard deviation which can be interpreted in terms of 
the time -dependent velocity U(t) as the root-mean-square value of (U(t) - U). The 
standard deviation is thus a measure of the magnitude of velocity fluctuations. The 
standard deviation, of course, gives no information on the temporal frequency spectra 
of the velocity fluctuations. 

The skewness of a histogram can be estimated by 



(Uj - Ua) 3 Cj 

a 2 


(B3) 


A histogram that is symmetric about the mean velocity has zero skew. A positive value 
of skew indicates that the histogram is weighted toward velocities higher than the mean 
or that it has a long tail above the mean. There is no familiar analogy to the histogram 
skew in the analysis of time -dependent functions. This is also true of the excess 
(or kurtosis) E. The excess is calculated by 


E = 


I 

I 

i=l 


(Uj - Ua) 4 Ci 


(B4) 


In the absence of velocity shear, the velocity in turbulent flow has often been char- 
acterized by the classical Gaussian, or normal, velocity probability density function. 

The normal probability density is symmetric so that the skew is zero. The excess is 
defined in such a way that a normal density has zero excess. Negative values of excess 
indicate a velocity probability density function that is more peaked or pointed than a nor- 
mal density function. That is, negative excess means that the velocities for the measured 
probability density, when compared to a normal distribution with the same standard devi- 
ation, are more closely bunched about the mean. Positive values of excess characterize 
a probability density that is flatter than a normal density function. 

The usual methods of analyzing time -varying signals (ref. 21), such as a power 
spectrum or autocorrelation, allow the identification of many common functions of time 
by characteristic shapes. For example, Gaussian white noise yields a power spectrum 
that is a constant (a straight line). The corresponding autocorrelation is a spike at a 
delay time of zero. A histogram of white noise also has a characteristic shape, that is, 
that of a Gaussian or normal probability density function. The interpretation of histo- 
grams is greatly aided by the ability to recognize some characteristic shapes and relate 
them to shapes of power spectra and autocorrelations. 
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Frequency Velocity 

Sketch (b) 


A velocity that varies sinusoidally with time, U(t) = U a + U r sin cut, is illustrated 
in sketch (b). The power spectrum is a spike at the frequency u>. The width of the spike 
is the frequency resolution A to, and the height of the spike is proportional to the value 
of U r . The shape of the autocorrelation is a cosine function of delay time. The delay 
t im e between the peaks of the autocorrelation is the period of U(t). The constant magni- 
tude of the cosine wave indicates that neither the frequency co nor the magnitude of U r 
varies with time. The shape of the corresponding histogram is symmetric about mean 
velocity U a . The width of the histogram is twice the velocity variation U r . The two- 
peaked (bimodal) shape indicates that the velocity is fluctuating between two values. 

The exact time variation that causes a measured histogram is difficult to determine 
for several reasons. Turbulence is present in many flow situations. For example, the 
square wave in sketch (c) has a low turbulence level superimposed on the high velocity, 
and a larger turbulence fluctuation exists at the lower velocity state. The resultant histo- 
gram resembles the sum of two normal probability density distributions. However, no 
reliable method of defining the two separate flow states by direct examination of the histo- 
gram has been found. 

The most useful method for the analysis of histograms of time -dependent velocities 
has been found to be an indirect method of three steps. The first step is to form an 
analytical hypothesis of the type of time variation that yields the histogram. The hypoth- 
esis accounts for as much of the time variation as possible by estimation of the physics 
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Sketch (c) 

of the flow situation. The remainder of the time variation may be accounted for by a 
Monte Carlo method. That is, one of the flow parameters is a random variable. The 
hypothetical model is then used to generate a time history of the velocity. 

The second step is to form an analytical model of the data-acquisition process so 
that a simulated histogram can be generated. The pertinent features of the data- 
acquisition and data-reduction process are modeled so that the simulated histogram will 
be identical to the measured histogram if there is no defect in the hypothetical flow model. 

The third step is to compare the simulated and measured histograms. If the com- 
parison is not satisfactory, the hypothetical flow model must be corrected. If the mea- 
sured and simulated histograms agree perfectly, the hypothetical flow model may not be 
accepted with complete assurance as representing the true time variation because other 
models could yield the same result. The fact that the time and frequency information is 
lost in the process of histogram construction precludes the possibility of eliminating all 
other time variations that may yield the same histogram. Therefore, care must be taken 
to assure that simulations satisfy the physical constraints of the flow. 
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TABLE L- LASER VELOCIMETER DATA AND FREE-STREAM FLOW CONDITIONS 


Run 

X, 

m 

y> 

m 

C, N, 

T r percent percent 
chord | chord 

“s> 

deg 

«w» u L,a» V L| a, U Ljm ,V LjmJI U L>1 ,V L}lf a u , ct v , 

m/sec m/sec m/sec m/sec m/sec m/sec m/sec m/sec 

S R,u 

s R,v 

E u 

E v 

D U 

D v 

VcOj 

m/sec 

Mco 

Rc 

Date 

1 

0.1362 

0.0981 

3 

4.01 

29.28 

0.00 

19.44 

114 

109 

125 

113 

102 

113 

24 

25 

-0.30 

0.26 

-0.36 

-0.40 

148 

141 

178.6 

0.493 

1.331 

Feb. 18 

2 

.1362 

.0874 

3 

4.02 

22.26 

.01 

19.44 

112 

103 

122 

113 

133 

93 

27 

34 

-.71 

-.41 

1.25 

.04 

169 

200 

178.6 

.493 

1.331 

Feb. 18 

3 

.1362 

.0824 

3 

4.02 

18.98 

.02 

19.44 

104 

104 

107 

103 

138 

93 

30 

25 

.18 

.82 

-.38 

.47 

289 

226 

178.6 

.493 

1.331 

Feb. 18 

4 

.1362 

.0773 

3 

4.02 | 

15.63 

.03 

19.44 

116 

104 

128 

111 

105 

87 

30 

31 

.33 

.22 

-.25 

-.34 

202 

187 

178.5 

.493 

1.327 

Feb. 18 

5 

.1364 

.0742 

3 

3.95 

13.60 

-.13 

19.50 

112 

131 

115 

165 

174 

175 

36 

40 

.56 

-.10 

-1.05 

-1.73 

419 

463 

170.6 

.490 

1.469 

Feb. 25 

6 

.1362 

.0723 

3 

4.02 

12.35 

.04 

19.44 

107 

117 

120 

124 

130 

165 

29 

34 

-.35 

.07 

-.32 

-.98 

288 

154 

178.5 

.493 

1.327 

Feb. 18 

7 

.1364 

.0694 

3 

3.94 

10.45 

-.15 

19.50 

145 

147 

174 

175 

184 

180 

41 

40 

-.56 

-.72 

-1.34 

-1.23 

359 

827 

170.6 

.490 

1.469 

Feb. 25 

8 

.1363 

.0673 

3 

. 3.98 

9.07 

-.08 

19.44 

104 

108 

102 

105 

130 

93 

28 

33 

' .70 

.47 

.48 

-.28 

200 

132 

178.5 

.493 

1.327 

Feb. 18 

9 

.1364 

.0643 

3 

3.94 

7.10 

-.16 

19.50 

55 

73 

74 

90 

81 

196 

35 

48 

1.20 

1.02 

1.92 

.53 

1747 

1438 

171.3 

.489 

1.449 

Feb. 25 

10 

.1364 

.0643 

3 

3.94 

7.10 

-.16 

19.50 

79 

83 

97 

93 

87 

87 

43 

38 

.15 

.11 

-.86 

-.17 

1668 

1778 

171.3 

.489 

1.449 

Feb. 25 

11 

.1364 

.0643 

3 

3.94 

7.10 

-.16 

19.50 

109 

132 

112 

170 

194 

196 

42 

48 

.72 

.11 

-.48 

-1.69 

470 

566. 

171.3 

.489 

1.449 

Feb. 25 

12 

.1363 

.0623 

3 

4.02 

5.88 

.01 

19.44 

84 

90 

97 

98 

97 

124 

34 

36 

-.16 

-.75 

-.31 

-.04 

133 

137 

178.5 

.493 

1.327 

Feb. 18 

1? 

.1364 

.0623 

3 

3.97 

5.88 

-.08 

19.50 

70 

103 

76 

186 

71 

198 

31 

68 

1.25 

.27 

3.89 

-1.46 

1333 

682 

171.2 

.490 

1.453 

Feb. 25 

14 

.1364 

.0603 

3 

3.96 

4.57 

-.08 

19.50 

64 

78 

76 

103 

71 

196 

40 

57 

1.44 

1.02 

2.54 

-.19 

842 

815 

171.0 

.490 

1.461 

Feb. 25 

15 

.1364 

.0583 

3 

3.96 

3.26 

-.09 

19.50 

67 

77 

79 

95 

71 

188 

42 

53 

1.43 

1.00 

1.88 

-.04 

1141 

1002 

171.0 

.490 

1.461 

Feb. 25 

16 

.1363 

.0571 

3 

4.02 

2.47 

.01 

19.44 

82 

77 

99 

95 

92 

103 

46 

38 

.30 

.15 

-.83 

-.61 

697 

524 

177.9 

.490 

1.316 

Feb. 18 

17 

.1364 

.0563 

3 

3.96 

1.95 

-.10 

19.50 

57 

67 

71 

82 

79 

196 

36 

45 

1.52 

1.34 

3.30 

1.47 

1398 

1285 

171.0 

.490 

1.461 

Feb. 25 

18 

.1364 

.0544 

3 

3.96 

.70 

-.11 

19.50 

51 

52 

74 

74 

79 

87 

32 

34 

.86 

1.08 

.64 

1.09 

1283 

1159 

170.6 

.490 

1.469 

Feb. 25 

19 

.1364 

.0534 

3 

3.95 

.04 

-.12 

19.50 

137 

156 

186 

191 

197 

196 

49 

45 

-.04 

-.71 

-1.62 

-1.02 

318 

288 

171.3 

.489 

1.449 , 

Feb. 25 

20 

.1363 

.0533 

3 

4.02 

-.02 

.01 

19.44 

93 

67 

140 

95 

145 

98 

58 

40 

.07 

.24 

-1.33 

-1.04 

365 

570 

177.9 

.490 

1.316 

Feb. 18 

21 

.1363 

.0533 

3 

4.02 

-.02 

.01 

19.44 

108 

84 

133 

98 

130 

93 

48 

38 

-.40 

-.18 

-.72 

-.67 

427 

489 

177.9 

.490 

1.314 

Feb. 18 

22 

.1225 

.0923 

5 

8.93 

26.32 

8.81 

19.50 

143 

110 

163 

108 

166 

108 

33 

13 

-.96 

1.16 

-.26 

6.60 

264 

366 

175.4 

.490 

1.420 

Feb. 24 

23 

.1235 

.0872 

5 

8.83 

22.92 

8.71 

19.50 

151 

110 

169 

111 

174 

108 

29 

13 

-1.02 

2.10 

-.08 

8.64 

231 

262 

175.4 

.490 

1.420 

Feb. 24 

24 

.1240 

.0823 

5 

8.95 

19.69 

8.83 

19.50 

144 

106 

169 

105 

174 

108 

34 

18 

-.62 

1.54 

-.99 

4.31 

395 

325 

175.4 

.490 

1.420 

Feb. 24 

25 

.1247 

.0774 

5 

8.97 

16.44 

8.86 

19.50 

124 

102 

133 

98 

174 

98 

33 

20 

.35 

1.93 

-1.17 

4.21 

364 

489 

175.4 

.490 

1.420 

Feb. 24 

26 

.1255 

.0723 

5 

8.97 

13.05 

8.86 

19.50 

125 

104 

138 

98 

176 

93 

36 

26 

.33 

1.47 

-1.31 

1.37 

768 

493 

170.1 

.491 

1.476 

Feb. 24 

27 

.1263 

.0675 

5 

8.94 

9.86 

8.82 

19.50 

125 

106 

145 

98 

184 

95 

43 

29 

.27 

1.38 

-1.39 

.93 

594 

289 

170.1 

.491 

1.476 

Feb. 24 

28 

.1262 

.0673 

5 

9.01 

9.74 

8.90 

19.50 

146 

95 

174 

93 

176 

98 

40 

22 

-.78 

2.24 

-.99 

6.42 

491 

384 

168.8 

.490 

1.507 

Feb. 25 

29 

.1270 

.0626 

5 

8.97 

6.61 

8.86 

19.50 

176 

107 

184 

98 

184 

93 

24 

35 

-2.79 

1.51 

7.17 

.91 

1738 

704 

168.8 

.490 

1.507 

Feb. 25 

30 

.1277 

.0602 

5 

8.78 

4.99 

8.64 

19.50 

184 

127 

189 

142 

189 

191 

21 

45 

-3.73 

.55 

14.55 

-1.44 

1870 

533 

168.8 

.490 

1.507 

Feb. 25 

31 

.1277 

.0592 

5 

8.87 

4.34 

8.74 

19.50 

170 

131 

194 

157 

194 

196 

44 

47 

-1.34 

.33 

.19 

-1.48 

1227 

679 

169.7 

.490 

1.487 

Feb. 25 

32 

.1282 

.0567 

5 

8.79 

2.67 

8.66 

19.50 

116 

100 

169 

118 

194 

103 

61 

50 

-.02 

.17 

-1.32 

-.61 

1188 

1368 

169.7 

.490 

1.487 

Feb. 25 

33 

.1283 

.0540 

5 

8.99 

.91 

8.88 

19.50 

31 

36 

38 

43 

24 

41 

17 

20 

1.19 

.94 

1.54 

.65 

683 

1413 

170.1 

.489 

1.474 

Feb. 25 

34 

.1091 

.0947 

7 

14.85 

29.60 

13.60 

19.50 

140 

122 

163 

124 

166 

124 

32 

12 

-.85 

-.24 

-.47 

7.34 

436 

590 

173.7 

.491 

1.431 

Feb. 24 

35 

.1102 

.0894 

7 

14.94 

26.05 

13.60 

19.45 

117 

109 

117 

105 

102 

103 

28 

21 

.45 

.79 

-.04 

1.27 

1004 

351 

179.3 

.491 

1.383 

Feb. 20 

36 

.1105 

.0848 

7 

15.36 

23.07 

13.85 

19.45 

101 

106 

102 

105 

102 

105 

27 

24 

.96 

.97 

.98 

1.83 

663 

176 

179.3 

.491 

1.383 

Feb. 20 

37 

.1124 

.0770 

7 

15.37 

17.80 

13.86 

19.45 

105 

112 

110 

111 

102 

103 

35 

30 

.49 

.82 

-.20 

.01 

807 

458 

179.0 

.490 

1.379 

Feb. 20 

38 

.1133 

.0751 

7 

15.14 

16.45 

13.72 

19.45 

114 

114 

117 

121 

97 

95 

34 

37 

.54 

.51 

-.65 

-.70 

1336. 

521 

179.5 

.491 

1.378 

Feb. 20 

39 

.1148 

.0673 

7 

15.37 

11.24 

13.85 

19.45 

134 

144 

156 

175 

189 

196 

42 

43 

.00 

-.14 

-1.37 

-1.50 

514 

819 

179.0 

.490 

1.379 

Feb. 20 

40 

.1157 

.0634 

7 

15.40 

8.62 

13.87 

19.45 

105 

107 

115 

118 

97 

108 

39 

38 

.20 

.12 

.09 

-.05 

913 

1157 

179.0 

.490 

1.379 

Feb. 20 


►P- 

tt- 


TABLE I.- Continued 


Run 

X, 

m 

y , 1 
m ! 

Tr 

c , 

percent 

chord 

1 N ‘ 1 
percent 

chord 

| 

deg 

deg 

u L,a> 

m/sec 

v L,a.|u L ,m.|vL,m,|uL,l> 

m/secl m/sec m/sec'm/sec 

V L,1»| a u> 1 a v» | 
m/sec|m/sec m/sec 

S R,u 1 

s R,v ! 

E u 1 

E v 

Du : 

Dy 

Voo, 

m/sec 

m to 1 

R c 1 

Date ! 

41 

0.1163 

0.0615 

7 

15.32 

7.31 

13.83 

19.45 

92 

98 

1 107 

108 

89 

98 

41 

36 

0.04 

-0.08 

-0.20 

0.17 

798 

1285 

179.0 

0.490 

1.379 

Feb. 20 

42 

.1174 

.0575 

7 

15.25 

4.59 

13.79 

19.45 

65 

71 

74 

82 

51 

51 

29 

33 

1.04 

.96 

1.54 

.79 

929 

1227 

179.5 

.490 

1.372 

Feb. 20 

43 

.1178 

.0555 

7 

15.21 

3.25 

13.82 

19.45 

51 

96 

57 

103 

52 

82 

21 

33 

.11 

.23 

-.39 

.58 

736 

811 

179.5 

.490 

1.372 

Feb. 20 

44 

.1186 

.0537 

7 

15.09 

1.98 

13.69 

19.45 

27 

36 

29 

47 

25 

51 

10 

22 

.32 

.92 

-.40 

.20 

723 

731 

179.5 

.490 

1.372 

Feb. 20 

45 

.1185 

.0518 

7 

15.44 

.79 

13.90 

19.45 

20 

22 

25 

27 

32 

28 

10 

11 

.57 

.52 

-.39 

-.50 

405 

380 

179.8 

.490 

1.363 

Feb. 20 

46 

.0881 

.0956 

9 

25.42 

33.98 

18.13 

19.44 

143 

113 

166 

113 

169 

105 

31 

29 

-.85 

.72 

-.13 

-.61 

115 

65 

177.0 

.494 

1.378 

Feb. 18 

47 

.0881 

.0956 

9 

25.42 

33.98 

18.13 

19.44 

105 

105 

107 

100 

94 

98 

33 

38 

.79 

.35 

-.34 

-.43 

153 

53 

177.0 

.494 

1.378 

Feb. 18 

48 

.0911 

.0866 

9 

25.37 

27.75 

18.11 

19.44 

145 

119 

169 

139 

176 

167 

38 

39 

-1.69 

-.39 

2.64 

-.23 

408 

302 

176.5 

.491 

1.361 

Feb. 18 

49 

.0942 

.0765 

9 

25.46 

20.81 

18.14 

19.44 

114 

115 

156 

162 

169 

173 

56 

52 

-.41 

-.20 

-1.06 

-1.05 

498 

662 

176.7 

.491 

1.354 

Feb. 18 

50 

.0974 

.0675 

9 

25.36 

14.56 

18.11 

19.44 

132 

123 

161 

136 

179 

131 

49 

44 

-.80 

-.65 

-.10 

.13 

381 

535 

176.7 

.491 

1.354 

Feb. 18 

51 

.1001 

.0575 

9 

25.70 

7.77 

18.22 

19.44 

81 

86 

99 

95 

102 

103 

35 

27 

.05 

-.05 

-.43 

1.39 

490 

797 

176.6 

.490 

1.348 

Feb. 18 

52 

.1016 

.0536 

9 

25.55 

5.03 

18.17 

19.44 

46 

51 

48 

54 

48 

43 

19 

22 

1.46 

.92 

4.60 

1.63 

303 

559 

176.6 

.490 

1.348 

Feb. 18 

53 

.1016 

.0536 

9 

25.55 

5.03 

18.17 

19.44 

45 

47 

51 

56 

57 

43 

18 

20 

-.10 

.10 

-.14 

-.31 

473 

671 

176.6 

.490 

1.348 

Feb. 18 

54 

.1029 

.0495 

9 

25.60 

2.21 

18.19 

19.44 

47 

46 

47 

55 

38 

46 

20 

22 

.83 

.39 

.03 

-.40 

447 

401 

176.8 

.490 

1.342 

Feb. 18 

55 

.1029 

.0495 

9 

25.60 

2.21 

18.19 

19.44 

29 

30 

32 

32 

33 

25 

13 

14 

1.33 

1.47 

2.66 

3. 36 

254 

415 

176.8 

.490 

1.342 

Feb. 18 

1 56 

.1029 

.0495 

9 

25.60 

2.21 

18.19 

19.44 

26 

27 

26 

28 

29 

24 

8 

8 

.38 

.45 

.69 

.50 

304 

437 

176.8 

.490 

1.342 

Feb. 18 

57 

.1029 

.0495 

9 

25.60 

2.21 

18.19 

19.44 

23 

23 

21 

23 

21 

19 

7 

5 

1.28 

1.29 

1.96 

1.77 

197 

280 

176.8 

.490 

1.338 

Feb. 18 

58 

.1064 

.0482 

9 

23.66 

.69 

17.54 

19.44 

14 

14 

18 

18 

21 

18 

7 

8 

.96 

.50 

1.28 

-.40 

34 

92 

176.8 

.490 

1.338 

Feb. IB 

59 

.1064 

.0482 

9 

23.66 

.69 

17.54 

19.44 

15 

11 

20 

16 

22 

24 

7 

7 

-1.02 

.19 

-.62 

-1.04 

172 

196 

177.3 

.490 

1.334 

Feb. 18 

60 

.0694 

.0935 

11 

35.78 

36.78 

20.79 

19.45 

107 

130 

120 

131 

163 

134 

40 

15 

.34 

-.36 

-1.10 

2.36 

1316 

479 

166.5 

.490 

1.529 

Feb. 23 ! 

61 

.0727 

.0845 

11 

35.85 

30.49 

20.80 

19.45 

147 

136 

156 

139 

158 

142 

26 

18 

-.86 

-1.32 

.12 

3.58 

519 

639 

166.6 

.490 

1.535 

Feb. 23 

62 

.0745 

.0798 

11 

35.84 

27.18 

20.80 

19.45 

126 

139 

140 

147 

153 

147 

33 

20 

-.01 

-.90 

-1.16 

.96 

1380 

1338 

166.6 

.490 

1.535 

Feb. 23 

63 

.0764 

.0752 

11 

35.75 

23.92 

20.78 

19.45 

142 

144 

153 

149 

153 

155 

28 

19 

-1.23 

-1.12 

.80 

1.56 

419 

378 

166.6 

.490 

1.535 

Feb. 23 

64 

.0780 

.0704 

11 

35.88 

20.60 

20.80 

19.45 

137 

137 

151 

149 

163 

155 

30 

29 

-.85 

-.96 

-.01 

.23 

661 

584 

166.8 

.491 

1.532 

Feb. 23 

65 

.0797 

.0657 

11 

35.92 

17.32 

20.81 

19.45 

131 

132 

143 

144 

153 

165 

29 

30 

-.41 

-.40 

-.57 

-.74 

823 

855 

166.8 

.491 

1.532 

Feb. 23 

66 

.0815 

.0611 

11 

35.89 

14.08 

20.81 

19.45 

95 

103 

112 

121 

97 

118 

40 

41 

-.14 

-.29 

-.44 

-.34 

1104 

917 

166.8 

.491 

1.532 

Feb. 23 

67 

.0832 

.0563 

11 

35.97 

10.74 

20.82 

19.45 

62 

59 

66 

67 

46 

46 

26 

24 

.94 

1.30 

1.34 

3.25 

1173 

773 

166.7 

.490 

1.529 

Feb. 23 

68 

.0841 

.0545 

11 

35.84 

9.43 

20.80 

19.45 

62 

61 

69 

62 

48 

46 

24 

28 

1.37 

1.27 

2.56 

2.36 

1616 

831 

166.7 

.490 

1.529 

Feb. 23 

69 

.0847 

.0527 

11 

35.89 

8.18 

20.81 

19.45 

50 

53 

56 

56 

46 

49 

20 

19 

-.03 

.25 

-.09 

.03 

1355 

828 

166.7 

.490 

1.529 

Feb. 23 

70 

.0858 

.0507 

11 

35.68 

6.70 

20.77 

19.45 

29 

30 

30 

32 

21 

23 

13 

12 

1.79 

1.41 

4.56 

3.31 

1081 

988 

166.9 

.490 

1.527 

Feb. 23 

71 

.0858 

.0507 

11 

35.68 

6.70 

20.77 

19.45 

24 

37 

30 

43 

24 

38 

12 

19 

.13 

.99 

-.73 

.69 

808 

838 

166.9 

.490 

1.527 

Feb. 23 

72 

.0869 

.0474 

11 

35.78 

4.42 

20.78 

19.45 

28 

29 

29 

31 

24 

23 

10 

11 

1.66 

1.80 

4.98 

5.51 

954 

1246 

166.9 

.490 

1.527 

Feb. 23 

73 

.0874 

.0451 

11 

36.00 

2.89 

20.83 

19.45 

17 

18 

22 

22 

19 

22 

10 

10 

.93 

.86 

.44 

.29 

665 

886 

166.7 

.490 

1.524 

Feb. 23 

74 

.0435 

.1024 

13 

47.89 

48.52 

22.79 

19.45 

129 

136 

153 

147 

163 

162 

33 

27 

-.49 

-.56 

-.91 

-.26 

569 

241 

180.3 

.491 

1.364 

Feb. 20 

75 

.0450 

.0988 

13 

47.90 

45.96 

22.79 

19.45 

116 

133 

13a 

144 

153 

155 

33 

27 

.04 

-.73 

-1.02 

-.16 

545 

234 

179.9 

.490 

1.360 

Feb. 20 

76 

.0467 

.0967 

13 

47.45 

44.26 

22.73 

19.45 

149 

138 

151 

142 

151 

147 

13 

18 

-2.63 

-1.50 

10.72 

3.34 

504 

212 

164.8 

.490 

1.561 

Feb. 23 

77 

.0467 

.0967 

13 

47.45 

44.26 

22.73 

19.45 

148 

138 

151 

144 

145 

147 

15 

19 

-2.29 

-2.45 

9.97 

6.82 

629 

596 

164.8 

.490 

1.561 

Feb. 23 

78 

.0466 

.0950 

13 

47.90 

43.26 

22.79 

19.45 

101 

126 

112 

144 

153 

155 

39 

36 

.39 

-.34 

-1.02 

-.82 

341 

240 

180.0 

.490 

1.361 

Feb. 20 

79 

.0479 

.0928 

13 

47.69 

41.60 

22.76 

19.45 

151 

135 

158 

147 

153 

155 

24 

28 

-1.54 

-.91 

2.86 

.23 

568 

259 

164.8 

.490 

1.561 

Feb. 23 

80 

.0525 

.0827 

13 

47.49 

34.32 

22.73 

19.45 

104 

104 

138 

142 

145 

147 

49 

52 

-.30 

-.44 

-1.06 

-1.15 

856 

616 

164.9 

.489 

1.551 

Feb. 23 

81 

.0544 

.0789 

13 

47.31 

31.54 

22.71 

19.45 

80 

84 

92 

95 

87 

87 

37 

38 

.59 

.53 

-.12 

-.36 

1385 

836 

164.9 

.489 

1.551 

Feb. 23 

82 

.0556 

.0751 

13 

47.54 

28.93 

22.74 

19.45 

52 

57 

61 

69 

48 

38 

30 

36 

1.45 

1.30 

2.19 

1.21 

1036 

874 

164.9 

.489 

1.551 

Feb. 23 

83 

.0573 

.0713 

13 

47.47 

26.20 

22.73 

19.45 

50 

59 

61 

77 

46 

46 

29 

37 

1.2? 

1.08 

1.68 

.67 

1121 ' 

921 

166.2 

.490 

1.537 

Feb. 23 

84 

.0587 

.0676 

13 

47.56 

23.61 ; 

22.74 

19.45 

29 

56 

39 

80 

19 

54 

20 

41 

1.31 

1.21 

1.24 

.81 

1198 

534 

166.2 

.490 

1.537 

Feb. 23 


TABLE I.- Continued 


Run 

X, 

m 

y. 

m 

C, N, 

T r percent percent 
chord chord 

“s. 

deg 

°w ^L,a» ^L,a> ^L,mj 
deg m/sec m/sec m/sec 

v L,m> u L,l’ V L,1> °u. 0 V’ 

m/sec,m/sec m/sec m/sec m/sec 

s R,u 

S R,V 

E u 

E v 

D e 

V 

D °°’ 

v m/sec 


R c 

Date 

85 

0.0603 

9.0639 

13 

47.53 

20.96 

22.74 

19.45 

30 

37 

37 

50 

19 

19 

17 

21 

1.46 

0.87 

2.01 

-0.03 

845 

549 

166.4 0.490 

1.535 

Feb. 23 

86 

.0603 

.0639 

13 

47.53 

20.96 

22.74 

19.45 

34 

48 

39 

60 

38 

62 

18 

17 

.95 

-.34 

.93 

-.43 

1106 

1074 

166.4 

.490 

1.535 

Feb. 23 

87 

.0619 

.0604 

13 

47.45 

18.44 

22.73 

19.45 

21 

39 

23 

53 

19 

56 

10 

18 

.56 

.59 

.45 

-.28 

1121 

1018 

166.4 

.490 

1.535 

Feb. 23 

88 

.0619 

.0604 

13 

47.45 

18.44 

22.73 

19.45 

25 

33 

25 

46 

19 

56 

12 

17 

1.99 

.76 

5.16 

-.16 

910 

723 

166.6 

.491 

1.536 

Feb. 23 

89 

.0628 

.0585 

13 

47.39 

17.06 

22.72 

19.45 

20 

35 

24 

50 

19 

58 

11 

17 

.46 

.80 

.03 

.01 

867 

724 

166.6 

.491 

1.536 

Feb. 23 

90 

.0650 

.0529 

13 

47.48 

13.11 

22.73 

19.45 

22 

37 

24 

51 

19 

56 

10 

18 

.54 

.53 

.09 

-.37 

1102 

682 

166.3 

.490 

1.532 

Feb. 23 

91 

.0675 

.0474 

13 

47.36 

9.15 

22.72 

19.45 

24 

56 

25 

62 

19 

62 

9 

11 

.54 

-.46 

.48 

.48 

886 

944 

166.3 

.490 

1.532 

Feb. 23 

92 

.0690 

.0437 

13 

47.39 

6.53 

22.72 

19.45 

23 

54 

23 

.62 

19 

62 

8 

13 

1.16 

-.03 

1.57 

.42 

669 

563 

166.3 

.490 

1.532 

Feb. 23 

93 

.0698 

.0419 

13 

47.36 

5.24 

22.72 

19.45 

22 

54 

22 

62 

19 

62 

8 

12 

.87 

-.29 

.00 

-.42 

661 

530 

166.2 

.489 

1.527 

Feb. 23 

94 

.0704 

.0401 

13 

47.45 

4.00 

22.73 

19.45 

23 

53 

21 

62 

19 

62 

6 

13 

1.66 

.18 

2.34 

-.27 

340 

418 

166.2 

.489 

1.527 

Feb. 23 

95 

.0714 

.0381 

13 

47.36 

2.53 

22.72 

19.45 

23 

51 

21 

60 

19 

62 

7 

13 

2.04 

.77 

3.94 

.90 

296 

475 

166.5 

.490 

1.529 

Feb. 23 

96 

.0721 

.0347 

13 

47.79 

.30 

22.78 

19.45 

22 

48 

21 

53 

19 

62 

6 

13 

1.65 

.38 

3.05 

.10 

154 

435 

166.5 

.490 

1.529 

Feb. 23 

97 

.0224 

.1018 

15 

59.55 

53.67 

24.10 

19.45 

147 

142 

151 

142 

153 

142 

18 

7 

-2.86 

-1.04 

9.80 

4.93 

1028 

572 

166.7 

.490 

1.524 

Feb. 23 

98 

.0265 

.0926 

15 

59.56 

47.07 

24.10 

19.45 

145 

143 

151 

147 

145 

147 

24 

15 

-1.67 

-2.19 

2.98 

8.46 

811 

463 

166.8 

.490 

1.522 

Feb. 23 

99 

.0307 

.0835 

15 

59.49 

40.49 

24.10 

19.45 

115 

137 

140 

147 

153 

147 

40 

23 

-.40 

-1.00 

-1.14 

.69 

1136 

958 

166.8 

.490 

1.522 

Feb. 23 

100 

.0348 

.0743 

15 

59.50 

33.88 

24.10 

19.45 

128 

136 

138 

144 

145 

155 

26 

26 

-.52 

-.67 

-.22 

.30 

951 

1137 

166.8 

.490 

1.522 

Feb. 23 

101 

.0393 

.0653 

15 

59.22 

27.28 

24.07 

19.45 

78 

88 

89 

103 

48 

87 

34 

42 

.86 

.72 

.13 

-.56 

1070 

1251 

166.8 

.490 

1.522 

Feb. 23 

102 

.0409 

.0608 

15 

59.46 

24.16 

24.09 

19.45 

68 

75 

79 

87 

51 

46 

32 

38 

1.04 

1.12 

1.05 

.71 

898 

1123 

166.8 

.490 

1.522 

Feb. 23 

103 

.0430 

.0561 

15 

59.46 • 

20.78 

24.09 

19.45 

35 

38 

38 

42 

21 

22 

18 

19 

1.34 

1.28 

1.75 

1.33 

1353 

1194 

166.9 

.490 

1.526 

Feb. 23 

104 

.0450 

.0517 

15 

59.44 

17.61 

24.09 

19.45 

39 

40 

46 

46 

38 

43 

20 

21 

.99 

.91 

.32 

.39 

1684 

1602 

166.9 

.490 

1.526 

Feb. 23 

105 

.0469 

.0471 

15 

59.53 

14.35 

24.10 

19.45 

43 

43 

49 

50 

43 

49 

21 

22 

.74 

.73 

-.13 

-.19 

821 

1260 

166.9 

.490 

1.526 

Feb. 23 

106 

.0490 

.0425 

15 

59.51 

11.03 

24.10 

19.45 

33 

34 

35 

37 

21 

24 

16 

16 

1.58 

1.52 

2.67 

2.66 

976 

1550 1 

166.9 

.490 

1.530 

Feb. 23 

107 

.0512 

.0379 

15 

59.42 

7.68 

24.09 

19.45 

38 

39 

42 

43 

28 

41 

17 

18 

1.12 

1.12 

1.18 

1.10 

920 

1609 

167.0 

.491 

1.536 

Feb. 23 

108 

.0519 

.0361 

15 

59.48 

6.42 

24.10 

19.45 

38 

40 

42 

45 

38 

49 

17 

18 

1.02 

1.01 

1.18 

.93 

961 

1658 

167.0 

.491 

1.536 

Feb. 23 

109 

.0529 

.0342 

15 

59.39 

5.01 

24.09 

19.45 

36 

37 

39 

42 

35 

43 

17 

17 

1.28 

1.15 

1.56 

1.40 

1101 

1586 

165.3 

.491 

1.558 

Feb. 23 

110. 

.0536 

.0325 

15 

59.43 

3.81 

24.09 

19.45 

33 

31 

34 

32 

24 

23 

17 

11 

1.82 

1.51 

3.34 

3.64 

1146 

1489 

166.4 

.491 

1.550 

Feb. 23 

Ill 

.0543 

.0307 

15 

59.49 

2.54 

24.10 

19.45 

46 

36 

62 

40 

81 

23 

24 

16 

.55 

1.22 

-.97 

1.66 

420 

325 

166.4 

.491 

1.550 

Feb. 23 

112 

-.0063 

.1030 

17 

74.98 

62.29 

25.38 

19.44 

119 

135 

143 

144 

153 

155 

38 

26 

-.19 

-.57 

-1.11 

.33 

484 

302 

176.6 

.492 

1.359 

Feb. 19 

113 

-.0020 

.0933 

17 

75.14 

55.33 

25.39 

19.44 

135 

132 

153 

147 

153 

155 

36 

32 

-.79 

-.53 

-.20 

-.49 

1053 

527 

176.6 

.492 

1.359 

Feb. 19 

114 

.0024 

.0843 

17 

75.07 

48.75 

25.39 

19.44 

130 

124 

138 

134 

153 

142 

27 

29 

-.36 

-.42 

.46 

-.04 

1181 

1177 

176.6 

.492 

1.359 

Feb. 19 

115 

.0066 

.0753 

17 

75.09 

42.23 

25.39 

19.44 

85 

89 

92 

95 

102 

98 

26 

29 

.44 

.73 

.25 

.29 

1136 

1239 

176.6 

.491 

1.353 

Feb. 19 

116 

.0066 

.0753 

17 

75.09 

42.23 

25.39 

19.44 

84 

88 

92 

95 

102 

98 

27 

30 

.27 

.60 

.35 

.33 

1139 

1242 

176.6 

.491 

1.353 

Feb. 19 

117 

.0086 

.0703 

17 

75.29 

38.69 

25.41 

19.44 

62 

63 

65 

67 

73 

93 

17 

18 

.27 

.00 

-.54 

-.22 

1350 

1335 

176.6 

.491 

1.353 

Feb. 19 

118 

.0107. 

.0662 

17 

75.19 

35.67 

25.40 

19.44 

45 

48 

51 

53 

48 

47 

18 

20 

.52 

.62 

-.18 

-.34 

608 

670 

176.3 

.490 

1.345 

Feb. 19 

119 

.0129 

.0623 

17 

75.02 

32.75 

25.38 

19.44 

35 

40 

38 

46 

25 

41 

15 

17 

1.21 

1.01 

1.69 

.67 

1094 

1516 

176.3 

.490 

1.345 

Feb. 19 

120 

.0151 

.0570 

17 

75.20 

28.99 

25.40 

19.44 

21 

22 

25 

25 

22 

24 

10 

10 

.32 

.61 

-.40 

.04 

721 

646 

176.3 

.490 

1.345 

Feb. 19 

121 

.0173 

.0532 

17 

74.97 

26.12 

25.38 

19.44 

20 

22 

23 

25 

19 

25 

9 

10 

.68 

.52 

.12 

-.20 

964 

959 

176.9 

.491 

1.342 

Feb. 19 

122 

.0204 

.0487 

17 

74.42 

22.58 

25.34 

19.44 

18 

23 

23 

25 

23 

25 

9 

9 

.77 

.57 

.25 

.03 

728 

873 

176.9 

.491 

1.342 

Feb. 19 

123 

.0214 

.0442 

17 

75.07 

19.63 

25.39 

19.44 

18 

22 

22 

25 

25 

25 

10 

9 

.72 

.58 

.01 

.07 

596 

814 

176.9 

.491 

1.342 

Feb. 19 

124 

.0237 

.0396 

17 

75.00 

16.25 

25.38 

19.44 

14 

24 

16 

27 

19 

25 

5 

10 

.10 

.47 

-.94 

-.22 

745 

825 

176.3 

.488 

1.334 

Feb. 19 

125 

.0258 

.0352 

17 

75.00 

13.06 

25.38 

19.44 

14 

25 

17 

27 

20 

25 

5 

10 

-.04 

.44 

-.98 

-.38 

382 

689 

176.3 

.488 

1.334 

Feb. 19 

126 

.0259 

.0352 

17 

74.94 

13.03 

25.39 

19.45 

18 

28 

19 

30 

19 

40 

4 

9 

-1.10 

.59 

1.12 

-.52 

505 

543 

169.5 

.491 

1.483 

Feb. 20 

127 

.0283 

.0297 

17 

75.06 

9.09 

25.39 

19.44 

14 

25 

16 

28 

20 

25 

5 

10 

.11 

.45 

-.95 

-.41 

506 

625 

176.3 

.488 

1.333 

Feb. 19 

128 

.0284 

.0298 

17 

74.97 

9.12 

25.39 

19.45 

13 

22 

17 

j 25 

20 

22 

5 

11 

-.05 

.64 

-.96 

-.36 

658 

760 

170.0 

.491 

1.475 

Feb. 20 

129 

.0310 

.0244 

17 

74.95 

5.19 

25.39 

19.45 

18 

23 

22 

27 

20 

19 

10 

11 

.82 

.51 

.16 

-.43 

407 

684 

178.1 

.490 

1.401 

Feb. 20 
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X, 

m 

y, 

m 
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c, 
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N, 
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deg 
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E u 

E v 

Du 

Dv 
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m/sec 


R c 

Date 

130 

0.0318 

0.0220 

17 

75.15 

3.54 

25.40 

19.45 

22 

25 

26 

30 

48 

25 

11 

12 

0.51 

0.29 

0.19 

0.73 

257 

597 

178.2 

0.490 

1.392 

Feb. 20 

131 

.0318 

.0220 

17 

75.15 

3.54 

25.40 

19.45 

29 

28 

32 

30 

19 

20 

10 

9 

.65 

.81 

-.84 

-.47 

283 

530 

178.6 

.490 

1.389 

Feb. 20 

132 

-.0444 

.0972 

19 

96.89 

69.91 

27.19 

19.50 

132 

143 

138 

147 

140 

147 

23 

15 

-1.11 

-1.91 

1.79 

5.42 

404 

766 

176.0 

.491 

1.411 

Feb. 24 

133 

-.0420 

.0928 

19 

96.81 

66.63 

27.19 

19.50 

130 

147 

138 

149 

133 

147 

25 

20 

-.59 

-1.18 

.27 

3.00 

508 

675 

176.0 

.491 

1.411 

Feb. 24 

134 

-.0396 

.0883 

19 

96.76 

63.28 

27.18 

19.50 

110 

143 

130 

155 

130 

165 

38 

31 

-.45 

-1.48 

-.27 

3.00 

588 

799 

176.0 

.491 

1.411 

Feb. 24 

135 

-.0375 

.0838 

19 

96.87 

60.02 

27.19 

19.50 

101 

138 

122 

152 

130 

165 

40 

33 

-.26 

-.92 

-.42 

1.02 

536 

649 

176.1 

.490 

1.404 

Feb. 24 

136 

-.0352 

.0793 

19 

96.87 

56.71 

27.19 

19.50 

91 

122 

97 

134 

94 

124 

32 

36 

.25 

-.37 

.42 

.22 

682 

788 

176.1 

.490 

1.404 

Feb. 24 

137 

-.0328 

.0750 

19 

96.77 

53.48 

27.18 

19.50 

70 

94 

76 

108 

71 

98 

29 

42 

1.07 

.47 

1.66 

-.55 

917 

1306 

176.4 

.490 

1.401 

Feb. 24 

138 

-.0305 

.0704 

19 

96.80 

50.10 

27.18 

19.50 

50 

87 

61 

126 

73 

175 

24 

52 

.50 

.41 

-.97 

-1.10 

662 

1884 

176.4 

.490 

1.401 

Feb. 24 

139 

-.0285 

.0661 

19 

96.91 

47.00 

27.19 

19.50 

63 

108 

75 

136 

87 

155 

22 

47 

-.18 

.01 

-.99 

-1.08 

681 

1802 

176.4 

.490 

1.401 

Feb. 24 

140 

-.0285 

.0661 

19 

96.91 

47.00 

27.19 

19.50 

92 

115 

115 

144 

97 

165 

43 

48 

.29 

-.24 

-.83 

-.97 

927 

1363 

174.3 

.490 

1.417 

Feb. 24 

141 

-.0261 

.0616 

19 

96.86 

43.65 

27.19 

19.50 

105 

124 

122 

149 

122 

180 

41 

47 

-.19 

-.58 

-.48 

-.50 

614 

1055 

174.3 

.490 

1.417 

Feb. 24 

142 

-.0239 

.0572 

19 

96.89 

40.42 

27.19 

19.50 

86 

no 

99 

134 

87 

175 

38 

47 

.57 

.01 

-.05 

-.93 

893 

1232 

174.3 

.490 

1.417 

Feb. 24 

143 

-.0216 

.0526 

19 

96.92 

37.05 

27.20 

19.50 

76 

101 

89 

118 

92 

95 

35 

45 

,B0 

.39 

.61 

-.79 

1220 

1165 

173.6 

.489 

1.420 

Feb. 24 

144 

-.0194 

.0482 

19 

96.95 

33.82 

27.20 

19.50 

37 

70 

48 

98 

23 

41 

23 

48 

.87 

.96 

-.10 

-.11 

1053 

1220 

173.6 

.489 

1.420 

Feb. 24 

145 

-.0169 

.0438 

19 

96.81 

30.50 

27.19 

19.50 

43 

81 

56 

98 

43 

41 

24 

47 

.46 

.76 

-.69 

-.37 

793 

902 

176.3 

.490 

1.402 

Feb. 24 

146 

-.0147 

.0395 

19 

96.81 

27.33 

27.19 

19.50 

41 

72 

55 

90 

41 

87 

24 

43 

.45 

.92 

-.74 

.14 

763 

901 

176.3 

.490 

1.402 

Feb. 24 

147 

-.0124 

.0348 

19 

96.87 

23.90 

27.19 

19.50 

37 

62 

44 

80 

43 

46 

19 

36 

.95 

1.16 

.53 

1.07 

856 

1309 

176.3 

.490 

'.402 

Feb. 24 

148 

-.0103 

.0304 

19 

96.96 

20.70 

27.20 

19.50 

39 

47 

44 

50 

41 

49 

17 

20 

.96 

.82 

.84 

-.05 

499 

783 

176.3 

.490 

1.402 

Feb. 24 

149 

-.0079 

.0261 

19 

96.85 

17.47 

27.19 

19.50 

36 

41 

39 

43 

47 

41 

16 

17 

1.20 

1.21 

1.88 

1.38 

514 

929 

176.3 

.490 

1.402 

Feb, 24 

150 

-.0057 

.0215 

19 

96.94 

14.13 

27.20 

19.50 

33 

38 

37 

42 

32 

43 

14 

16 

1.24 

1.33 

2.21 

1.80 

775 

1264 

176.3 

.490 

1.402 

Feb. 24 

151 

-.0034 

.0172 

19 

96.88 

10.93 

27.19 

19.50 

31 

33 

33 

36 

23 

25 

12 

13 

1.40 

1.34 

2.86 

2.36 

550 

988 

176.3 

.490 

1.402 

Feb. 24 

152 

-.0016 

.0138 

19 

96.85 

8.40 

27.19 

19.50 

21 

22 

24 

24 

20 

24 

9 

10 

.66 

.76 

-.13 

-.06 

373 

719 

170.1 

.491 

1.479 

Feb. 24 

153 

.0002 

.0099 

19 

96.97 

5.59 

27.20 

19.50 

18 

19 

20 

22 

16 

20 

9 

9 

1.19 

1.01 

1.02 

.72 

196 

416 

170.1 

.491 

1.479 

Feb. 24 

154 

-.0566 

.0964 

T.E. 

113.97 

72.01 

19.50 

19.50 

114 

151 

122 

160 

130 

165 

29 

21 

.28 

-1.30 

-.52 

1.77 

873 

1142 

170.8 

.490 

1.469 

Feb. 24 

155 

-.0535 

.0869 

T.E. 

114.13 

65.45 

19.50 

19.50 

113 

146 

117 

160 

97 

165 

28 

30 

.47 

-.65 

-.23 

-.61 

797 

538 

170.1 

.490 

1.483 

Feb. 24 

156 

-.0500 

.0776 

T.E. 

113.99 

58.94 

19.50 

19.50 

89 

112 

102 

149 

97 

165 

38 

49 

.31 

.02 

-.29 

-1.37 

326 

962 

170.1 

.490 

1.483 

Feb. 24 

157 

-.0466 

.0681 

T.E. 

113.96 

52.31 

19.50 

19.50 

87 

108 

89 

134 

87 

170 

23 

45 

1.09 

.30 

2.76 

-1.23 

1845 

847 

169.0 

.490 

1.509 

Feb. 24 

158 

-.0432 

.0587 

T.E. 

113.91 

45.76 

19.50 

19.50 

85 

115 

87 

147 

87 

170 

23 

46 

.78 

.06 

2.96 

-1.32 

2215 

679 

169.0 

.490 

1.509 

Feb. 24 

159 

-.0403 

.0491 

T.E. 

114.21 

39.18 

19.50 

19.50 

65 

72 

76 

90 

76 

46 

31 

42 

.48 

.99 

.69 

.22 

1047 

1572 

169.0 

.490 

1.509 

Feb. 24 

160 

-.0365 

.0398 

T.E. 

113.89 

32.60 

19.45 

19.45 

49 

60 

57 

65 

44 

46 

21 

17 

.50 

.25 

-.56 

-.83 

903 

1974 

164.7 

.490 

1.572 

Feb. 23 

161 

-.0365 

.0398 

T.E. 

113.89 

32.60 

19.45 

19.45 

43 

47 

48 

53 

39 

41 

19 

21 

.84 

.76 

.04 

-.33 

1519 

1588 

164.7 

.490 

1.572 

Feb. 23 

162 

-.0334 

.0304 

T.E. 

114.03 

26.10 

19.45 

19.45 

40 

45 

44 

47 

38 

43 

17 

19 

.96 

.88 

.54 

.13 

1206 

1331 

164.7 

.490 

1.572 

Feb. 23 

163 

-.0316 

.0257 

T.E. 

113.94 

22.80 

19.45 

19.45 

40 

44 

44 

47 

48 

43 

16 

19 

.93 

1.03 

.58 

.53 

889 

1057 

164.9 

.491 

1.571 

Feb. 23 

164 

-.0299 

.0210 

T.E. 

113.91 

19.52 

19.45 

19.45 

36 

42 

39 

46 

38 

43 

15 
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Figure 2,~ Two-component fringes in sample volume. 
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Figure 3.- Laser and optics table installation in plenum chamber of Langley high-speed 7- by 10-foot tunnel. 
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L-76-1348 

Figure 7.- Model wing mounted in Langley high-speed 7- by 10-foot tunnel looking downstream. 
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Figure 9.- Locations of pressure transducers and tufts. ID denotes identification. 
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Figure 15.- Velocities upstream of stagnation point of hemisphere at 0° angle of attack, 0.55 Mach number, 

and Reynolds number of 1.8 x 10®. 





Figure 17.- Diagram of component fringe pattern in sample volume. 
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Figure 20.- Histograms above transducer 5. Concluded. 
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Figure 23.- Histograms above transducer 11. 
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Figure 23.- Histograms above 
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Figure 24.- Histograms above transducer 13. 
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Figure 26.- Histograms above transducer 17. Continued. 
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Figure 26.- Histograms above transducer 17. Concluded. 
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Figure 30. - Histogram moments for scan above transducer 5. 
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Figure 32.- Histogram moments for scan above transducer 9. 
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Figure 34.- Histogram moments for scan above transducer 13. 
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Figure 36.- Histogram 
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(a) Measured delay times between various pairs of transducers. 
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(c) Frequency of pressure pulse passage found from autocorrelations. 
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(d) Spacing of pressure pulses from ratio of Vp to J2. 

Figure 46. - Correlation summary from online measurement of first positive 

correlation peak. 
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Figure 47. - Model of vortex stream used to generate simulated histograms. 
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Figure 48.- Measured (left 2 columns) and simulated (right 2 columns) histograms 

above transducer 7. 
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Figure 49. - Measured (left 2 columns) and simulated (right 2 columns) histograms 

above transducer 5. 
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Figure 52.- Correlations and Fourier transforms of correla 
tions for pressure at transducer 5. Concluded. 
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Figure 54.- Correlations and Fourier transforms of 
correlations for pressure at transducer 9. 
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Figure 55. - Correlations and Fourier transforms of correla 
tions for pressure at transducer 11. Concluded. 
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Figure 56.- Correlations and Fourier transforms of correla 
tions for pressure at transducer 15. Concluded. 
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Figure 57.- Correlations and Fourier transforms of 
correlations for pressure at transducers 20 and 2. 
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(g) Transducer 15. 
Figure 58.- Continued. 
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Figure 60.- Map of two components of standard deviation of velocity corrected for velocity bias 

and scaled by tunnel velocity. 







(a) Theoretical histogram of Gaussian distribution. 



(b) Effect of measuring absolute velocity values above 5 m/sec. 
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(c) Additional effect of high-pass filter. 

Figure 64.- Effect of data processing on mean velocity and standard deviation. 
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(b) Values calculated from histograms generated from profile 1. 

Figure 65.- Comparison of measured histogram moments to moments simulated from hypothetical profile 1. 
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(Ab) Simulated histograms 
with Ul = Vl = 42 m/sec 
and ct u = a v = 18 m/sec. 



(Ac) Simulated histograms 
with Ul = Vl = 42 m/sec 
and cr u = a v = 18 m/sec. 
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and cr u = u v = 22 m/sec. 
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and o u = cr v = 22 m/sec. 




(Ca) Measured histograms at 
y = 0.0532 m. 


0 150 0 150 

V^, m/sec m/sec 




(Da) Measured histograms at 
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(Cb) Simulated histograms 
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and ct u = a v = 22 m/sec. 
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(Db) Simulated histograms 
with Ul = Vl = 5 m/sec 
and cr u = a v = 22 m/sec. 
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(Dc) Simulated histograms 
with Ul = Vl = -5 m/sec 
and cr u = ct v = 22 m/sec. 


Figure 66.- Measured and simulated histograms above transducer 17. 



percent 


Measured Profile 1 Profile 2 



(Ea) Measured histograms at (Eb) Simulated histograms (Ec) Simulated histograms 

y = 0.0442 m. with Ul = Vl = 1 m/sec with Ul = Vl = -12 m/sec 

and ct u = cr v = 22 m/sec. and ct u = cr v = 22 m/sec. 



(Fa) Measured histograms at (Fb) Simulated histograms (Fc) Simulated histograms 

y = 0.0396 m. with Ul = Vl = -1 m/sec with Ul = Vl = -18 m/sec 

and ct u = ct v = 22 m/sec. and a u = cr v = 22 m/sec. 



(Ga) Measured histograms at (Gb) Simulated histograms (Gc) Simulated histograms 

y = 0.0352 m. with Ul = Vl = -3 m/sec with Ul = Vl = -19 m/sec 

and a u = cr v = 22 m/sec. and ff u = cr v = 22 m/sec. 



V^, m/sec U L , m/sec V^, m/sec U^, m/sec V^, m/sec U^, m/s ec 


(Ha) Measured histograms at (Hb) Simulated histograms (He) Simulated histograms 

y = 0.0352 m. with Ul = Vl = -3 m/sec with Ul = Vl = -19 m/sec 

and cr u = cr v = 22 m/sec. and cr u = cr v = 22 m/sec. 


cn 

cn 


Figure 66.- Continued. 



percent 


Measured 


Profile 1 


Profile 2 


cn 


03 



(la) Measured histograms at 
y = 0.0298 m. 



(Ja) Measured histograms at 
y = 0.0244 m. 



(Ka) Measured histograms at 
y = 0.0220 m. 
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(La) Measured histograms at 
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(Lb) Simulated histograms 
with U L = V L = -1 m/sec 
and a u = a y = 22 m/sec. 



(Ic) Simulated histograms 
with Ul = Vl = -17 m/sec 
and a u = cr v = 22 m/sec. 
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and cr u = cr v = 22 m/sec. 
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(Lc) Simulated histograms 
with Ul = Vl = -8 m/sec 
and o u = ct v = 22 m/sec. 


Figure 66.- Concluded. 
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(a) Comparison of values calculated from measured histograms at transducer 17 with 

arbitrarily chosen profile 2. 
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(b) Values calculated from histograms generated from profile 2. 

Figure 67. - Comparison of measured histogram moments to moments simulated from hypothetical profile 2. 
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